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ABSTRACT. Many problems in probability connected with m events occurring randomly 
in space can be substantially simplified by dealing first with an indefinite number of events 
distributed in a Poisson series. A combinatorial problem arising in cosmic-ray counting 
discussed recently by Schrédinger provides a good example, and several properties of the 
corresponding distribution are derived. A generalization of this problem is considered, and 
the method is also applied to the random division of a line. 


$1. INTRODUCTION 


N a recent note Schrodinger (1951) has considered a combinatorial problem 
] arising in connection with cosmic-ray counters when the possibility of more 

than one event occurring in each counter is allowed for. This is one of a 
class of problems in which it is substantially simpler to consider the number of 
events occurring as indefinite, and distributed in a Poisson series. The solution 
when the number of events is fixed can then readily be derived as a corollary, 
and properties such as the mean and mean square deviation can be simply deduced. 


§2. THE COSMIC-RAY PROBLEM 


The problem put to Schrédinger was as follows: 

I. A tray of r equal counters is exposed to a cosmic-ray shower and hit by 
n(21) rays. A particular ray hits one, and only one, counter, and the probability 
of its hitting any particular counter is 1/r. What is the probability of precisely 
m counters going off? 

We shall consider the following alternative problem: 

II. Cosmic rays occur randomly in space and time, the mean number passing 
through a tray of r equal counters in a given time being A. Each ray hits one, 
and only one, counter. What is the probability of precisely m counters going off ? 

Denote the probability in I by w(n; 7, m), and in II by W(A; 7, m). The 
solution of I can easily be derived from II. It is well known that in II the 
probability of precisely m rays passing through the counters is e *A"/n!. Hence, 

PROC. PHYS. SOC. LXV, 5—A 21 


306 C. Domb 


we can divide W(A; r, m) into mutually exclusive cases in which precisely rays 
have passed through, and write identically 


W(A; 7, m) = = e aN (ns 7th). ~~“ eoeae (1) 


Thus, if we expand W(A; r, m)e’ as a power series in A we can obtain w(n; 1, m). 
The solution of II is very straightforward since the counters are independent. 
The probability of no ray passing through a particular counter is equal to 
exp (—A/r)(=q). The probability of one or more rays passing through is 
[1—exp(—A/r)](=p). The probability of m counters going off is 


! 
(Ce pg cS = om) eee (2) 


m\(r—my)! 
and the distribution is binomial. 

The practical problem, to estimate A or p from m, is a straightforward problem 
of binomial estimation. As pointed out by Schrédinger, there is a difficulty in 
deriving an estimate when m=r. It is clear that in this case, if we set some 
arbitrary probability limit, the only conclusion we can reasonably draw is that 
G<Qo OrA>Ap. In statistical terminology we would say that the upper confidence 
limit has become infinite. 

To derive Schrédinger’s solution from (2) we must expand 


Cn (1—e7 Uym(e— Ary et ="C (ee —1)™@ eae (3) 


as a power series in A, and pick out the coefficient of fie !, This is easily seen to be 


2 m re (Mt) r| 2 (m—t)” 

Con 2 EY face Tears: ae ‘Goi (4) 
which is equivalent to eqn. (6) of Schrodinger’s paper. It is clear from (3) that 
the series (4) is identically zero if n<m. 7 

The moments and moment generating functions can be derived in an 
analogous manner. Denoting, for example, the moment generating functions 
for the distributions I and II by f(n; 7, 0) and F(A; 7, 6) respectively, we have 


f(n;7, 0)= X wn; 7, mye 

Tae Pie Sines Cae Cae OTD A SE (5) 
and FA;7, = 2D” WA;r, m)e™. 

m=0 

eg ay An 
Hence, from (1) F(A; 7, 0) =" Se = ny FO On ae eee (6) 

n=0 . 
and a corresponding result holds for any of the moments. In_ particular, 
since for [ = 5% 

m=rp=r(1—e-*"), 

and im =rp? + rpg =P (1 —e- 4a pea ent) i geome o (7) 


we deduce that, for II 
m =r[1—(1~1/r)"], 
and m=r9[1 —2(1 —1/r)"+(1—2/r)"] +rf(1—1/r)"—(1 2p, fF ©) 


Needless to say it would have been a much more complicated matter to have 
derived (8) directly from (4). ; 


2+ = ew 9 ae 
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$3. THE DISTRIBUTION OF CLUSTERS 
A natural generalization which suggests itself is to divide the counters into 
groups through which exactly k rays have passed. Let us denote the probability 
that there are m counters of this type by w,(n; r, m) for I, and W,(A; 7, m) for II. 
Clearly 


fo 8) 
w(n; 1, m)= & w,(n; 1, m), 
k=1 


WA;7,m)= X W,(A; 17, m). 
k=1 


The solution of II is binomial as before, except that now the probability of 
k rays passing through a particular counter is exp (—A/r)(A/r)"/k!(=p,). Hence, 


WAST, m) = "Copa (Gp le py) tee (10) 
Expanding e*W,,(A; 7, m) as before we find that w,(n; 7, m) is the coefficient 
of A” /n! in 
A dkm ar Ae rem 
Cn ee a yee E ¥ al reat C » 
which, after a little simplification, gives 
Braise ee te ol n! 5 
ee Goh ck ee a) 


(where x >km and 7 is the integral part of n/k—m). 
Again the mean and mean square of the W,,(A; 7, m) distribution are given by 


= r Xr E —Alr 

na 9 aoe * 
ee (iNew oe A OTA oe (ie A ae a 
eter fig so kb \r) RUNG i | 


Proceeding as before, we deduce the corresponding moments of the w,(; r, m) 
distribution, 


mares) (= aes ee eee | (14) 
2 tatty "(2)" +m) O-)| 


The task of deriving these moments directly would have been even more: 
formidable than before. 

A number of authors have considered the problem of the distribution of 
aggregates in a random set of points (Silberstein 1945, Berg 1945, Mack 1948). 
The definition of an aggregate was as follows: a region of given size and shape. 
is moved all over the total length or area, and any set of k points falling within 
the region constitutes a k-aggregate. It is clear on this definition that any point 
may belong to several different aggregates, and for some purposes this may not’ 
be satisfactory. The present problem suggests an alternative possibility of 
dividing the total length or area arbitrarily into 7 cells, the number of points 
falling into each cell constituting an aggregate. ‘The distribution of aggregates 
is then given by eqns. (12) and (14) above. 
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§4. THE RANDOM DIVISION OF A LINE 


Suppose n points are taken at random on a line of length J. They subdivide 
the line into (n+1) intervals x,, x)... . #,4,;, and problems connected with the 
distribution of intervals have been discussed by many writers (see Moran 1947 
for references). ‘The method described in the present paper was first used for 
problems in which each point is the left-hand end-point of an interval of fixed 
length (Domb 1946). It can readily be applied to several other problems of a 
similar type. We always consider the associated problems in which events occur 
randomly on the line the mean number per unit length being A; the probability 
of m events occurring on the line is then exp(—A/)(Al)"/n!. Let us first apply 
the method to the determination of the probability distribution of the length 
of the first interval x, For the associated problem we have 


W(A; x) meer Neg ee (15) 
ee Nee 
Hence, expanding W(A; x,) e” as before, we find 
n(l—x,)"-1 n x,\%1 n>1 
w(n; x,) dx, = a a( ~ *) dx, (ey) Seok (16) 


It is clear that the second line of (15) can be ignored if we are interested in 7 >1. 
Similarly, if we require the joint probability distribution of the first 7 intervals, 
we have 


WONG Kis Xp 2 ON eae 


Sr eM ens rn) daa, ave Ti(a Ee ee ee ees (17) 
and hence 
LONE EN he eX) 
n(n—1)....(n—r+1) eat. eA ee 
= Sie [ae See (1 3 ered) (n 27). eee ewe (18) 
The moment generating function F(A; 0,) of the W(A; x,) distribution is 
: A 
U0, — A) —2%,(A— 0, = Gr Ne to 
‘ + re 2 de gh et cree eee (19) 
Hence, using an expansion analogous to (6), 
: q n'\ 10, 120," ge Pa 
f(s) = aga| ¢ lam aia ee aot |: tre (20) 


and in particular 
se! a! s_ irr! 
+ (a+)? “Tt  (n41)(n+2)? “1 (241)... (n tr)’ 


We can proceed similarly to determine the correlations between successive 
intervals. For example, in the Poisson problem, for a line of indefinite length 
x, and x, are uncorrelated. ‘The finite length of the line only introduces a 
correlation when less than two events occur. Hence, if we ignore this latter 
possibility and proceed as before, our results will be correct if m>2. For the 
Poisson distribution with a line of indefinite length we have the very simple result 


ris! 


The Use of a Random Parameter in Combinatorial Problems 309 


Thus, for the case of 2 points we deduce that 
— - | 
X4Xy = ——-.— , 
we (n+1)(n+2) 
ore ris! ere Fish itrs 
(r+s)! (n+r)...(n+7r+s) 
Prien berse a Wy Veet el [ese 


Seay er he aa Se ae kar,) 


xy"... xt= 
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ABSTRACT. A set of nuclear parameters (ranges and depths of potential wells) yielding 
the established values of the binding energy, quadrupole moment and magnetic moment 
of the deuteron, and consistent with the other experimental data relevant to the ground 
state of the neutron—proton system has been found. 


§1. INTRODUCTION 
HE ground state of the neutron—proton system is of particular interest in 
| the study of nuclear forces for the following reasons: it is a state of a 
two-particle system, so the analytical difficulties characteristic of many-body 
problems do not arise, and the complication of many-body forces is also avoided ; 
it is a stable nuclear state, so several different types of experimental data are 
available; and its several properties formerly led investigators to incompatible 
estimates of the range of the neutron—proton force. 
The main data relevant to the neutron—proton interaction are: 
(1) The binding energy E of the deuteron (E=2-23 Mey). 
(2) The quadrupole moment Q of the deuteron (O=2-76 x 10-?’ cm’). 
(3) The magnetic moment » of the deuteron (4 =0-8565 n.m.). 
‘4) Neutron—proton scattering cross sections and angular distributions. 
(5) Cross section for photodisintegration, and cross section for the converse 
process of deuteron formation by proton—neutron collision. 
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Of these, the first two can be determined quite accurately from experiment, 
and are also quite satisfactorily derivable from theory. The experimental 
accuracy of the third datum is also good, but theoretical calculation is made 
uncertain by relativistic effects ; however, it provides an estimate of the proportion 
of D state to be expected in the deuteron ground state—usually taken as 
4 pet -cent. 

In neutron—proton scattering, states other than the ground state are involved, 
and this means that assumptions about exchange properties which are not 
involved for (1), (2) and (3) have to be made to interpret data (4). In the case 
of low-energy scattering, an explicit exchange assumption can be replaced by the 
assumption of the charge-independence hypothesis, viz. that the 'S neutron— 
proton force is equal to the 1S proton—proton force as derived from proton—proton 
scattering experiments. Data of type (5) of useful accuracy have only recently 
become available (Snell, Barker and Sternberg 1949, Wilson, Collie and Halban 
1949, Barnes, Stafford and Wilkinson 1950). Datum (2) demands the inclusion 
of tensor forces in the interaction, but the quadrupole moment is not so large that 
the forces can be entirely tensor in character. 

Therefore the simplest type of investigation pretending to some degree of 
completeness is an examination of the ground-state of the neutron—proton 
system, central and tensor forces being considered: the postulated forces must 
lead to the values of FE, Q and pu given by experiment, and must be consistent 
with data of types (4) and (5). 

It is well known that the exact law of radial dependence of the nuclear 

interaction has little effect on most of the measurable quantities, at least at low 
energies. ‘This fact is implicit in the now well-established use of ‘effective 
range’. It is, therefore, quite sensible to make calculations using potential 
wells of any convenient shape, say square wells, find the appropriate ranges and 
depths, and translate into ‘effective ranges’ or into the equivalent parameters 
for other well-shapes as required (cf. Blatt 1948). 
_.. For the ground state of the deuteron, using a chosen well-shape and central 
and tensor forces, the nuclear parameters may be taken as four in number—two 
ranges and two depths—1.e. one parameter more than can be determined by using 
data (1)-(3). Rarita and Schwinger (1941) found for a square-well potential 
and equal tensor and central ranges a unique solution satisfying (1)-(3). The 
value of the common range was 2:8 x 10°'3cm. This solution also gave extremely 
good agreement with the data of type (4), and, with some reservations, those of 
type (5) then available. More recent experiments, particularly those on the 
scattering of neutrons in ortho and para hydrogen, however, have indicated that 
this range of 2:8 x 10cm is much too large to be compatible with data (4). 
Blatt and Jackson (1949) gave 1:56 x 10-'3cm (corresponding to a square well 
of roughly 1-7 10cm) as the best value for the effective range in the triplet 
state. ‘This figure has been revised to 1:74 x 10-cm by Hughes, Burgy and 
Ringo (1950). 

The results also of Barnes, Stafford and Wilkinson (1950) have been shown 
by Schiff (1950) to be consistent with an effective range of about 1-74 x 107% cm 
for the central force. 

Interesting figures have been given by ‘T’si-Ming Hu and Massey (1949) who 
derive optimum ranges to fit separately (a) the magnetic moment, (6) 1-6 Mev n—p 
scattering, (¢) scattering of thermal neutrons by hydrogen molecules, 
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(d) diffraction of neutrons by hydrogen in crystals, (e) photodisintegration of the 
deuteron. For a square well only the Rarita~Schwinger figure (for (a)) lies 
outside the region 1-4 to 2 x 10-8 cm. 

On the other hand, Bohm and Richman (1947) quote unpublished work by 
Sen (supported by later evidence) showing that 2:4 x 10-!%cm is the smallest 
range for which quadrupole moment and magnetic moment of the deuteron 
can be simultaneously satisfied. Both of these quantities are primarily dependent 
on the tensor force. 

The above mutually exclusive estimates of range are based on the assumption 
of a single range common to all the types of force which may be acting. The 
discrepancy of the results may therefore be read to mean that two or more ranges 
must be used if the properties of the neutron—proton ground state are to be 
satisfactorily accounted for. In the present work this hypothesis is investigated. 
We shall discuss the possibility of describing the ground-state phenomena in 
terms of a square potential well for the central force and a square well, but of 
different depth and width, for the tensor force. 

To conform to the indications of scattering and _ photodisintegration 
experiments we shall make the central force range of the order of magnitude 
1-4 to 2x 10-'%cm, while the tensor range and the two depths will be chosen 
to fit the data (1), (2) and (3). This is the main object of the following work. The 
general mode of interdependence of parameters and physical quantities will 
also be illustrated, and the usefulness and limitations of the differential analyser 
as a tool for this type of problem will be seen. 


32. THE DEUTERON EQUATIONS 
The equations for the neutron—proton system can be written 


a —«2)u= — B,w 


dw 6w 
EAP Sop (4,-0- 5) = —Bwu 
where 4,, x? and B,/2,/2 are proportional to the ponte force potential, binding 
energy and tensor force potential respectively. 
We shall use the square-well potentials 
A.=A, rea A,=05 ea 
BB rao Beard: 
In this paper the units of A, B, a, b are those which were found convenient in 
the course of calculation. ‘These units are defined, and the notation put into 
conformity with previous work, by the relations 
a=0-1157) where 7)=range of ordinary forces in 10-3 cm 
b=0-115¢ 7) where er) =range of tensor forces in 10° cm 
A=1-:83J where J =depth of ordinary well in Mev 
B=5:17yJ where yJ =depth of tensor well in Mev. 
The boundary conditions to be satisfied by equations I are: u, w both zero at 
r=0, u, w both +0 at r+, u, w, du/dr, dw/dr continuous at r=a, r=. 
Given a, 6, A, B, the two simultaneous differential equations for u, w with 
two-point boundary conditions (0, 00) for each of u, w would lead to an eigenvalue 
problem in «?. Since «” is in fact known from experiment (binding energy) 
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we shall instead treat equations I as an eigenvalue problem in one of the other 
parameters a, b, A or B, choosing whichever is mathematically more convenient. 

It is clear that the functions u, w, have different analytical forms in each of 
the three regions: (1) 0<r<a, (2) a<r<b, (3) b<r<, and that in region (1) 
a series solution will be tolerably convenient, while in region (3), where 
A,=B,=0, the equations separate and must have the solutions 


u=Uexp(-—ar), w= W(1+ 5 ar 2 ae x) exp (—ar7), 
where U, W are constants. 

To use a series or other analytic form of solution in region (2) and to find 
the analytic conditions for simultaneous fitting at r =a, b looks unwieldy, although 
Demeur (1950) succeeded in doing some work along these lines, viz. for ry =2°8 
and «=0-8, 0-9 and 1. This work enabled him to show that for 7)7=2-8 the 
experimental values of binding energy, quadrupole moment and magnetic 
moment could be accounted for only if «=1. 

In the present work fitting at y=a, 6 was obtained by linear combination of 
trial solutions. In region (2) numerical solutions were used, and, except for the 
final results, in which maximum accuracy was desired, the Cambridge differential 
analyser was used to perform this numerical integration. 

To make good use of this mechanical aid it is necessary, first, to obtain a 
fairly good estimate of the eigenvalue before attempting the exact solution and 
to have a rough idea of the values of uw, w and their derivatives in the numerical 
integration region, and secondly to treat A, B as given parameters (since they must 
be represented by gears on the machine) and either a or 6 as eigenvariable 
(preferably the greater, 5), since a, 6 can easily be varied by running the machine 
for a longer or shorter time before altering the gears representing A,, B,. 

The procedure used was as follows: 

(i) Approximate eigenvalues were found by an ordinary variational method. 
To be more explicit, a, b, A were allowed to take values of the sort of magnitude 
to be expected, 1.e. 8-7a=ry~1-2 to 3 (10° cm), 8-7b =er9~2-0 to 7 (10- cm), 
0:5484 =J~0 to 100 (Mev) and for each set of a, 6, A the corresponding 
eigenvalue B was found. ‘This approximate calculation also provided a rough 
guide to the behaviour of u, w, du/dr, dw/dr, and, consequently, an order of 
magnitude guide to the values of quadrupole moment and percentage D state 
resulting from each chosen set of a, b, A. 

(it) The exact solution could now be tackled. Sets of a, A, B were chosen 
on the guidance of the variational calculation; 6 was yet to be determined exactly. 
The differential analyser was set up for solving equations I with these values 
of A, B. It was also arranged that the value of A, could quickly be altered from 
A to 0 at the appropriate point (r =a) of the run. 

(iii) Since the differential analyser cannot be used at or too close to the 
singularity at the origin, a series solution was used for small values of r (7 <0-1 
was chosen for simplicity). The differential analyser started from the series 
solution values of u, w, du/dr, dw/dr at r=0-1, and the mechanical solution was 
run out to r=a, here A, was put equal to zero and the solution continued to 
the region r~b. 

(iv) The exact value of b was determined as described below, the corresponding 
eigenvalues u, w, du/dr, dw/dr at r<b deduced, and the asymptotic solutions 
fitted on to these differential analyser functions at r=). 
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§3. DETERMINATION OF 6 

The general solution of the two coupled second-order differential equations I 
involves four arbitrary constants, but the boundary condition u(0)=0, 
w(0)=0 reduces this number to two, and the homogeneity of the equations 
means that a third is an irrelevant multiplicative factor. Hence if (wu, #,), (U2, ®2) 
are two independent solutions non-singular at the origin, the eigensolution will 
be of the form (uw, + Buy, w, + Bws). 

With £6 chosen so that there was rough agreement with the variational 
estimates of the functions and their derivatives, a solution was run out to the 
region r~b. ‘The value of r at the point where the function u would fit on to 
the required asymptotic form U e~” we call b,, i.e. we define b, by 


dey 
Udr } 5a ree 


Similarly we define 6, from the differential analyser solution for w by the 
requirement 


EDp Sole Bay bear 
w dr}, — ‘ ay ott oes 0 8 ee 3a 


For an eigensolution 8 must be adjusted until 6, =b, (=6). At the expense 
of one or two extra runs on the machine (which were, however, very useful checks 
on mechanical and operational accuracy), solutions so close to the eigensolution 
were obtained that b could be determined as the point of intersection of graphs 
of 5, against 6 and 4, against f. 

Some comments will now be made on the several parts of the procedure 
outlined above. 


§4. THE VARIATIONAL CALCULATIONS 

Since the object of the variational calculation was rather to obtain information 
over a wide field than to get accurate results in a limited one, the simple polynomial 
forms u=r+¢'r? + f’r?, w=y'r? +8’r? were taken, these promising a reasonable 
reward for the amount of labour involved. (The primes distinguish the 
coefficients from «’s and f’s used elsewhere.) 

The above forms were fitted on to 

ar} _ 3 3 
u=Ue™,w=We™ (14 —+ =n) 
Ore. oer 
at r=b(>a) (or in one case r=a>5). 

This requirement completely determines the shape of the w function, only 
its size relative to u being adjustable. However, as experiment requires only 
about 4°% D state, the above forms, giving the greater flexibility to the wu function, 
seem reasonable. 

The information about eigenvalue behaviour obtained from this approxima- 
tion is shown in figs. 1 (a), (5). 

The circled points are accurately known, those on the B=0 axis from simple 
calculation, others from data given by Rarita and Schwinger. ‘The full curves 
are those deduced from the variational calculation; the dotted curves are 
estimated correction curves adjusted to pass through the accurately known 
points. It will be seen that the variational calculation looks tolerably good for 
b of the same order as a, but bad for 6 much greater than a. 
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Figure 2 shows the quadrupole moment behaviour indicated by the 


variational calculations. It is to be expected that these curves are much less 


reliable than those of fig. 1. 


J (Mev) 


b=02, 03217 
—— Values as calculated 
----Corrected values 
—:—~Comparison curves 

or 6=05,07,09 


YT (Mev) 


0 10 20 30 40 50 60 


b= 05,07 " 
—— Values as calculated 
---Corrected curves 


(0) 


{J (Mev) 


Figs. 1 (a) and(b). Ranges and well depths giving correct energy (rough variational calculation). 


In figs. 1-7: 
a=0-115 7) where ro>=range of ordinary forces in 10-18 cm. 
b=0-115 ev) where €/)= range of tensor forces in 10-13 cm. 
A=1-826 J where J=depth of ordinary well in Mev. 
B=5:17 yJ where y.J =depth of tensor well in Mev. 
Q=electric quadrupole moment in 10-2? cm?. 
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§5. THE MECHANICAL SET-UP ON THE DIFFERENTIAL ANALYSER 
In designing a set-up for the differential analyser the following considerations 
are important: (i) The functions on the displacement tables should be such 
that it is easy to satisfy the initial and boundary conditions. (ii) The generation 
of awkward multiples and awkward functions (e.g. (A,—«? — B,//2)w and 6w/r?) 
must be allowed for. (iii) Each integrator should achieve its maximum displace- 
ment at some point of the run. (This requires a knowledge of the magnitudes 
J (Mev) 
50 60 


4 =03217,05 


‘ 
—— Values as calculated 
—-- Corrected curves 


Q (10°“’cm?) 


Fig. 2. Quadrupole moments corresponding to ranges and depths from figs. 1 (a) and (8). 


of u, w, etc., conveniently provided by the variational calculation already carried 
out.) Also the multiples used on each shaft must be such that all rotate at 
similar speeds. ‘These conditions are required for efficient running and best 
accuracy. ‘The first of the above requirements was satisfied by using the 
integrated form of the equations 


du 4 of? S “7 
= +(4,-2) | udr=—B, | wdr, 


dw ["6w Pd a i 
=-| ag r+ (4,-« _ 5) wdr=—B, | udr, 


The integrator displacements involved were then u, w, du/dr, dw/dr, initial 

conditions for which were easily available, and the continuity conditions for 

these functions at r= a, ry =) meant that no recalculation or resetting was necessary. 
The integration, then, followed the scheme: 


. du 
es Aiea ) | udr *, 
ii aes ' dw 

B, | udr ——, 

dr 

{ ¢ 6w dw 

|e a 

Hale «oli A gpa | war eee 
dr rich 4/2 ; dr’ 
B, | war a 
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The essential integrations du/dr+u—> Ju dr, dw/dr>w- | w dr used four of the 
eight integrators of the Cambridge machine. 

The term {(w/r?)dr demanded the further use of one integrator if an input 
curve were used (say an input curve of 1/r, so that [(w/r?)dr= — Jwd(1/r)), or 
of two if it were desired to avoid the use of an input table. The second of these 
alternatives was chosen, and was carried out by means of a feedback set-up 
based on the two relations 


[-Kya([ Zar) = | war-K [ ar. Pes: (1) 
fo-Kya([ Sar) = [Zax [Sar Er (2) 


Thus, in step (1), given Jwdr and a displacement (r—K), we generated 
{(w/r) dr, and, similarly, from (2), having J(w/r)dr and displacement (r—K) 
we found | (w/r?) dr. 

This left two integrators free out of the eight: these were conveniently used 
for the multiplications by (4,—«?) and (A, —«?—B/1/2) which otherwise could 
only have been represented by the nearest fraction obtainable from a chain of 
gears of the ratios supplied for the differential analyser. 

Since no integrator was left for the B multiplications, the values of B were 
chosen so as to be easily obtainable from a chain of gears. 

The set-up, conforming to the requirements of (iii) above, and allowing 
for a running time of about half an hour per solution, was designed by 
Mr. B. Noble, to whom my thanks are due for much help and advice. 


§6. RESULTS OBTAINED FROM THE DIFFERENTIAL 
ANALYSER SOLUTIONS 

The results obtained from the differential analyser consist of the eigenvalue 
deduced for 5 given a, A, B, and tables of values of u, w, du/dr, dw/dr, for 
0-1 <r <b, the interval of tabulation being 0-025. The initial and asymptotic 
forms of the functions can then be calculated, and the corresponding quadrupole 
moment Q and the percentage D state D can be found. 

These, the most direct data obtainable, are shown (expressed in terms of 
To J etc.) an tablest, 

Our main interest, however, is not to find Q and D for arbitrary a, A, B, 
but to investigate the possibility of satisfying the simultaneous requirements 
E=2-23 Mev, O=2-76 x 10-?7cm?, D~4°% for a chosen value of a, the values 
of A, B, b being deduced. (At the time when this work was done the value of E 
used was 2-18 Mev and that of Q 2-73 x 10-?7cm®. The calculations have been 
revised to the Bell and Elliot value of H =2-237 mev, but further revision to the 
latest values of & and Q seems unnecessary. ‘The (small) corrections: necessary 
for adjustment to any new E, Q values can easily be estimated from the graphs 
given in this paper.) 

‘To obtain the sort of results mentioned above it was necessary to interpolate 
from a scattered lot of only about twenty eigenvalue sets, with three independent 
variables. However, the sets seemed capable of yielding very plausible results 
by means of the graph systems shown in figs. 3-7. These graphs probably 


become progressively less reliable as points on later graphs are deduced from 
the earlier ones. 
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The process works as follows: 

(i) A ‘preliminary eigenvalue plot’ (fig. 3) is made of B against 6. ‘This 
consists of a number of approximately parallel curves on each of which a, A 
are constant. From these the values of A for constant b can be read off, and a 
“full eigenvalue plot’ (fig. 4) of exactly the same type as was made for the 
variational calculation can be produced. The advantage of this procedure is 
that it allows us in our later graphs to reduce the number of variables; but, in 
contrast to the situation which would have arisen had the actual runs been 
performed with more constancy in the choice of parameters, we retain highest 
accuracy in the regions in which we are most interested. 


Table 1 
ro (10-cm) J (Mev) yJ (Mev) er9(10-8cem) Q(x10-*’cm?) D(%) 
1-31 46-4 14:5 Dit 2:64 5:07 
dod 46-4 15-4, 2:64 2°53 4-76 
1-74 27:3 sis DIST — — 
1-74 27-3 14:5 2:62 — = 
1-74 28-4 12-5; 2:84 2-76 4-79 
1-74 28-4 1315 2:69 2-52 4-46 
1-74 29°5 12-5, BET) 2°54 4-14 
1-74 29°5 113305) 2:63 2°36 4-19 
1-74 26°2 Sie5 2°81 2S) 5:16 
2:18 20°5 11-6 BOUT 2:62 4-29 
2-18 20:5 135 2°50 2:18 3-69 
2-61 15:8 10-6 DTS) 2:63 3-84 
2°61 15-8 11-6 2-64 2°41 B50 
2°8 14-0 10-8 ZA9) 2-70 3-89 
1:74 42-6 2:8: 5:30. SO7Al 2:64 
1-74 42-6 Sila 4-70 3-18 DeBXe 
2-18 26-2 Sl, 6-80 6°71 6-20 
2:18 29-5 Silis 3-92 2:92 2-04 
2°8 18-8, 2285 4-88 10-91 2:87 
2°8 18-8, 3-1, 4°37 8-04 2:03 
2°8 18-8, SEDs 3-87 9-63 Da5 2 


(ii) For the quadrupole moment a ‘preliminary Q plot’ of Q against 5 
analogous to the preliminary eigenvalue plot (i.e. 4, a constant on each curve) 
was made. The parameter B is, of course, uniquely determined for each point 
on these O curves by the full eigenvalue plot. Then a ‘full Q plot’ (fig. 5) of 
Q against A, each curve referring to constant a, b, is derived. From this a curve 
complying with Q =2-73 can be drawn for each a. ‘The ‘Q compliance curve’ 
(fig. 6) gives the relationship between 6 and A necessary for correct values of 
E and Q when B is taken at the value given by the eigenvalue plot. For each 
a also, a ‘D compliance curve’ for 4°%, D state can be drawn (fig. 6). Thus at 
the point of intersection of the two compliance curves for any particular a we can 
read off the 6, A, and hence deduce the B which will simultaneously give 
correct E, Q and D for the deuteron. 

It will be noted that the final value of B cannot be found directly from fig. 4, 
since 6 will not in general be one of the values chosen there. But an auxiliary 
curve of the same type and corresponding to the required value of 5 can readily 
be drawn for the purpose of this final determination. 
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Figs. 3 and 4. Ranges and well depths giving correct energy (differential analyser calculation). 
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Fig. 7. Variation of tensor range with ordinary range, the E, Q, D conditions all being satisfied. 
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The variation of 6 with a when all the deuteron conditions are satisfied is 
perhaps of interest. ‘This is shown in fig. 7. The Rarita-~Schwinger point 
here corresponds to D =3-9; the Guindon points to D~ 4-44, while the full curve 
is for D=4. In figs. 3-7, values obtained by Rarita and Schwinger, by Guindon 
(1948) and by Demeur (1950) have been included wherever comparison is 
possible. ‘The good agreement of these points gives grounds for a degree of 
confidence in the above method in spite of its dependence on extrapolation and 
graphical methods. 

At this stage the error was assessed at about 5°. A second set of differential 
analyser runs was done with values of a, b, A, B concentrated in the regions of 
interest. Interpolation to Q=2-73 x 10-?’cm? was performed with the aid of 
figs. 3-7, and table 2 obtained. 


Table 2. Final Parameter Sets deduced from Differential Analyser Solutions 


7) (10-18 cm) 1-4 1:6 1-8 2-0 

er) (10-48 cm) 2-94 3-04 3-06 3-04 
J (Mev) 46-86 39-33 32:33 26-19 
yJ (Mev) 11:57 10-04 9-52 9-52 
D(%) 4-62 4-14 4-00 4-00 


§7. ACCURACY AND FURTHER NUMERICAL WORK 

Since the precise amount of the error of the figures quoted above was very 
much in doubt—many of the tests carried out during the running of the machine 
were far from reassuring—numerical integrations using a desk calculating 
machine were carried out for six cases, viz. for 7) = 1-4 x 10-18 (two values of yJ/), 
1-6 x 10-8 (two values of yJ), 1-8 x 10-! and 2-0 x 10-% cm. 

The series solution was used to r=ry; step by step integration was performed 
with intervals of 0:1 x 10-13cm in 7, and the usual checking of fourth differences 
applied. The eigenvalues were found from two independent solutions, and the 
quadrupole moment and °% D state calculated using the optimum 5-, 6- and 
8-point formulae given by Bickley (1939). ‘The results were as shown in table 3. 


Table 3. Direct Results of the Computed Solutions 


r) (10-19 cm) 1-4 1-6 1-8 2-0 1-4 1-6 
€r,(10-2* cm) 3-121 3-108 3-005 3-010 2-692 2-968 
J (Mev) 51-78 39-60 32-23 26-19 51-78 39-60 
yJ (Mev) 9-52 9-52 9-52 9-52 11-42 9-96 
Q(10-27 cm?) 2-728 2-781 2-605 2-687 2-095 2-567 
D(%) 4-043 4-144 3-850 3-966 3-384 3-979 


It is of interest to see what errors this implies in the differential analyser 
results. The largest errors indicated in the differential analyser solutions are 
of the order of $°%, for ery and of the order of 1% for J, 79, yJ being presumed 
given. An error of 4% in yry may give rise to errors of 13% in Q, D, while an 
error of 1°, inJ may make Q, D in error by nearly 4%. Naturally, a small error 
in J may correspond to a large error in the binding energy—for r)=0-14 an 
error of 1%, in J is 0-5 Mev which is approximately 25°, of the binding energy, 
and although this is not an indication of the probable error in E£, it shows that an 
error of several per cent would not be surprising. 

The above error estimates apply to the differential analyser solutions.. The 
errors in the numerical solutions are very much less. Comparison with the 
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solution by the same method of the hydrogen atom ground-state equation 
d®y/dr? =x —2y/r, whose exact solution is y=re’, shows that the rounding-off 
errors are not likely to start building up in the deuteron equations in the region 
r<yr, with which we have been concerned. The results of table 3 may be 
reliable to 0-1%. 

Once again using figs. 3-7 to correct values to Q=2:-73 x 10-*” cm?, we finally 
obtain table 4. 


Table 4. Final Parameter Sets satisfying E =2:237 mev and Q =2:-73 x 10-?? cm? 


79 (10-18 cm) 1-4 1-6 1-8 2-0 
er) (10-1 cm) 3-12 3-08 3-06 3-04 
J (mev) 51-78 39-79 31-94 26-12 
yJ (MeV) 9-520 9-520 9-520 9-520 
D(%) 4-04 4-06 4-00 4:00 


It would appear safe to say that for E=2-237 Mev, O=2-73 x 10-?’ cm? this 
set of values is reliable to $°% in each of «79, J where 79, yJ have been chosen 
exactly ; and the corresponding values of D are not in error by as much as 4%. 

As mentioned before, the more recent values of the binding energy and 
quadrupole moment (Newell 1950) do not necessitate material alteration of the 
results given above. It is found that for E=2-22 Mev the quantities <7), J, D in 
table 4 are all changed by less than 0-2%%, i.e. by less than the suggested error 
magnitude. It therefore seems unnecessary, at least until values are well 
established, to give adjusted figures. Table 4 can therefore be taken as 
embodying the main results of this investigation. 
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ABSTRACT. An expression is given for the potential differences set up in colloidal 
suspensions of solid particles by the passage of sound waves which is correct to the first 
power of the zeta potential of the particles. "The question was examined in detail by Enderby 
in 1951. For a certain range of values of relevant parameters it was found necessary to intro- 
duce the simplification that the mobilities of the various ionic species in the electrolyte 
surrounding the colloidal particles were identical. In the present paper it is shown how the 
potential may be calculated quite generally without this restriction by the use of vector 
methods. 


§1. INTRODUCTION 

N a recent paper (Enderby 1951, to be referred to as E) a theoretical analysis 
] was given of the electrical effects accompanying the propagation of sound 
waves through suspensions of charged colloidal particles. It was shown in 
this paper that differences of electrical potential arise between different points 
in the wave train owing to the disturbance of the ionic double layers surrounding 
the particles; the spherical symmetry of the double layer field is modified by 
the flow of the liquid. Solutions for the potential distribution in the suspension, 
correct to the first power of the zeta potential of the particles, were obtained. 
These solutions are exact provided that the quantities (wk7T/c)!* are much 
greater or much less than the double layer thickness. The w, denote the ionic 
mobilities, o/27 is the frequency of the sound waves, and 7 the temperature. 
In the intermediate region of values it was found necessary to introduce the 
simplification that all the ionic mobilities were equal. ‘The following note 
indicates how the general solution can be evaluated for this region even though 
the ionic mobilities are not necessarily all identical. Hence the problem may 
now be regarded as solved over the whole range of frequency and double layer 
thickness. The method uses the concept of the ‘ion space’ and algebraic methods 
similar to those employed by Onsager and Fuoss (1932) in their treatment of 

the conductivity, viscosity and diffusion of mixtures of strong electrolytes. 


§2. GENERAL SOLUTION FOR THE POTENTIAL 


The method described in E was to evaluate the alterations in the mean 
potential and mean ionic concentrations in the neighbourhood of a single 
particle due to the passage of the sound waves. From these quantities the 
average potential in the wave train can easily be found by integration. 'To solve 
the equations of motion and of continuity for the ions it is necessary to expand 
the changes in potential and concentration in powers of the total charge on each 
particle. Thus, if ¢ denotes the alteration in the mean potential and p, the 
alteration in the concentration of the ion of type 7, at a point distant r from the 
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centre of the particle, the radius vector making angle @ with the direction of 
propagation of the sound waves, we write 


dur Ae Sa SO eB Nene g y=) a) ee (2.1) 
»v=1 n=0 


PHA Zea Oars Ar) aees 0), Ce een. (2.2) 
v=1 n=0 


where A is the amplitude of the sound waves, Qe the total charge on the particle, 
and P,, (cos @) is the usual Legendre polynomial. TheA, ,,and 7, , , are functions 
of r only. For a solution correct to the first power of Q or of the zeta potential 
we only need to consider those A and z functions for which the suffix » equals 
unity. All of these for which the second suffix is not equal to unity vanish, 
for reasons given in E. Hence the basic equations to be solved for a result 
correct to the first power of the zeta potential are eqns. (3.9) and (3.11) of E. 
When the ionic mobilities are not all equal, it is not in general possible to 
combine the two equations into a single fourth-order differential equation for 
the potential coefficient A,,, as in eqn. (3.12) of E. Suppose, however, we 
associate s potentials A, , , with the charge densities 7, , ; by means of the relations 


Nese yal ee oe eae alee (24) 


s denotes the number of ionic species present, « the dielectric constant of the 
electrolyte and A, is the operator defined in eqn. (3.10) of E. We now find, 
instead of the single eqn. (3.12) of E, the set of s equations: 


10 4ne'n;z; & 
Ay | (4: ere) ee kT ee B5Aj,3, | 


47e7n,2;, ay: ees 

= era | Vee +7, div v’ |. Sia aperak (2.4) 
In these equations we have used the notation of paper E. mn, denotes the 
concentration of the ion of type z outside the double layer region of the particle, 
and 2, its valency. V,, I’, and V’ are defined in eqns. (2.2), (3.9) and (2.3) 
respectively of E. Equations (2.4) give, when integrated, 


io 4n7e*n,z,; $ 
( — 5Te,) haa RT 2023 


J 


47e7n ;2 2 1 Dave 
~ (RT )%an(1 +5) | RG, x) + 100, y)+ (=) H(b, + | seid sed (2.5) 


All the quantities R(b, x), G(b, y), H(b, y) and g have been calculated in E. The 
boundary conditions used in eqn. (2.5) are, first, that X;,1,1 tends to zero at infinity, 
and secondly, that the field inside the particle is not affected by the conditions 
outside (eqn. (4.5) of E). 

To solve eqn. (2.5) we now introduce a space of s dimensions. The 
s potentials Ay 14, Ax 4,1--.-Ag4,1 may now be regarded as defining a vector I in 
this space. Similiarly let z have components 2,, 2,....%, and let t have 
components 1121/w1, Ne%/wWy,.... Equation (2.5) may now be written 


(eee RTM | RG, Aue 2160, ae (=): H(b, th Ae sQD18) 


Sound Waves in Suspensions 423 


where 4 is the matrix operator with components 


a 108 i 47e2n 12 i335 
aj kTw; eRT . @oesee 


The real parts of the characteristic numbers of the matrix A are all positive. 
This can easily be shown by first symmetrizing the matrix and then examining 
the quadratic form corresponding to the new matrix. Let B denote the 
transform of A defined by 


Boe SAS Paes | 2 piper ges... (2.8) 
where the components of the matrix S are given by 
Se Olas, ee Me ee Aa (2.9) 


We now find that if b,; is a component of B, 
; at 106 ;; 47e?(nn,) 12,2; 
aj ae k To; ek OF Ce On a ems Yer Wt a 


Hence the matrix B is symmetric; it is easily verified that the quadratic form 
corresponding to the real part of B is positive definite. 

Having obtained this result on the characteristic numbers, it is now possible 
to solve eqn. (2.6) for the potential vector by Onsager and Fuoss’s method. 
The solution will consist of a part decaying exponentially with increase of r, 
and a contribution varying as 7. The second part only is of interest for our 
purpose, and its value is 


bE li eee re 2 
(cRT«)(1 +b) 


From the definition of the A,,, it follows that the potential A, , is simply the 
scalar product of the two vectors z and 1. Hence, the effective dipole moment 
of the field round the particle is easily deduced from eqn. (2.11). The calculation 
of the potential difference ® across the half wavelength in the sound wave now 
follows exactly as in the earlier paper. The final expression for ® is 


___ 8NAAQ(1~p) 


ie AxLttat nen (2.11) 


(z. A“t)(G + (2r/ed)2H) 
RT(1+5) % ma2] 1+ 2p toy toy* | 


It only remains to evaluate the scalar products (z. A“!t) ; this can in principle 
be done directly, but the inversion of the matrix A is rather tedious. The 
procedure of Onsager and Fuoss could be applied, but the most convenient 
technique seems to be the application of the Lagrange—Sylvester expansion 
for A. This method was used by Guryanova (1936) to simplify the 
computations of Onsager and Fuoss. We have 


S (A—x,1)...(A—x,_,1)(A— 4,431) ...(A—%,1) 


A t= ap RE 2A 
pai Xp (%p— Hy) 00s (Hp —%y-1) (pm Hy t1).-- (%_—%s) ( ) 
where 1 denotes the unit matrix and x, x.,....%, are the characteristic 

numbers of the matrix A, that is, roots of the equation 
det (A= xia Ole ine = ge Mees (2.14) 


Except for the case of a binary or ternary electrolyte, a numerical solution 
of eqn. (2.14) is necessary. 


22-2 


324 F. Booth and F. A. Enderby 


§3. POTENTIAL SOLUTION FOR A BINARY ELECTROLYTE 


For a binary electrolyte s equals 2, and so an explicit expression for ® may 
be given without much difficulty, since the roots of the characteristic equation 
are easily determined. ‘The roots are 


1 al ca Yaa eel Lia ‘| 
smn 3[¢+ ap (G+ an) * + erp Las ~ a] atom 


é ar (= é =) P| 244, baal Sate (3.1) 


Using eqn. (2.13), we find 


DIG a ae ee ee (3.2) 
et —1¢)\ — to(m 421? + M229") 3 3) 
and for (z, A t) (z34 t) — RT xX Oy os) BU let ieetis ( i 


This reduces to the result already obtained when w, is set equal to a. 

The real part of the scalar product is very complicated, but gives a simple 
formula if z, and z, are equal in magnitude. A straightforward calculation 
shows that for the case 

BS Bon) gay Coe Ae eee (3.4) 
eqn. (2.12) may be written in the form of eqn. (3.27) of E multiplied by 


a factor B*, where 
% 2 fast | ct le Ee (3.5) 
=k" / Ww Reese SET} | at scl aaie : 


— 1/1 1 
and @® =43(w, +), 1 SS (- + ) oe DW ie euee taseire (3.6) 


Thus, in this simple case the more exact theory merely changes the f of E, 
defined in eqn. (4.7), to 8* of eqn. (3.5) above. Some idea of the extent of the 
modification can be obtained by comparing the two quantities 8 and 8*. If we 
attempt to use the approximate theory of E, §4, for the case s equal to 2 but 
, and w, unequal, by identifying the w of eqn. (4.7) with the mean mobility a, 
we And BRS en ae (3.7) 
This ratio is always less than unity; this means that the correct dotted curve 
of E, fig. 2, should lie below the one drawn there. For water the ratio B*/B is 
0-93. The maximum difference between 8* and B would be expected for an 
electrolyte in which one ion has a very high mobility and the other a very low one. 
An example of this type is acetic acid, which gives the value 0-36 for B*/B at 25°c. 

The general conclusion would seem to be, therefore, that the more exact 
theory gives potentials rather less than the approximate theory, but the effect is 
not important unless the mobilities differ considerably. In the latter case, 
however, it would be essential to use the more exact calculation in a careful 
comparison of experiment and theory. It is of interest to note that the previous 
calculations gave results rather greater than the experimental values, but for 
the reasons given in E it is not possible to conclude whether the present work 
agrees substantially better with the data. 
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Discontinuities in the Variation of Magnetization with 
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ABSTRACT. It is demonstrated that magnetic discontinuities accompany changes of 
temperature at constant field. The effect is investigated for hard-drawn iron and nickel and 
is accounted for in terms of two mechanisms, one associated with the decrease of coercivity 
with rise of temperature and the other associated with the disperse field due to the thermal 
oscillations of the carriers of magnetic moment. 


SNL RO DUE TR LON 
T has recently been shown experimentally by Tebble and Wood (unpublished 
i work) that changes of magnetization with temperature at constant field are 
reversible when the specimen is on a reversible portion of its hysteresis loop. 
When the constant field is in the region of the coercive force, however, 
irreversible changes of magnetization occur. It would be expected that the 
irreversible magnetization changes due to a change of temperature would take 
place in discontinuous jumps, as in the ordinary Barkhausen effect. Such 
discontinuities, which del Nunzio (1930) attempted but failed to detect, have 
now been demonstrated, and have been investigated in a qualitative manner for 
two materials. 
§2. EXPERIMENTAL 
The specimens, in the form of wires of 26 s.w.g. and of 17 cm length, were 
placed in a narrow glass tube, surrounded by a search coil. The temperature 
could be raised by the passage of hot air through the tube, and measured by 
means of a thermocouple placed in contact with the specimen. Any sharp 
discontinuous changes of magnetization were picked up by the search coil 
which was connected to a detecting system similar to that used by Tebble, 
Skidmore and Corner (1950). ‘The search coil and specimen were surrounded 
by a polarizing coil, and the whole system was magnetically screened. The 
minimum discontinuities detectable by the apparatus are of the order of 
1074 e201 
The first specimen, of hard drawn nickel, had a coercivity of about eight 
oersteds. No effect was detectable for this specimen in a saturating field of 
eighty oersteds, even with a temperature increase of as much as 100°c from 
room temperature. ‘This was also the case at remanence. When the specimen, 
after having been brought into a cyclic condition, was subjected to a demagnetizing 
field of nine oersteds, however, a rise in temperature of anything over 0:2°c 
produced many thousands of magnetic discontinuities, in such a direction as to 
increase the magnetization. When the temperature after having been raised 
was allowed to fall a very few discontinuities, still in the direction of increasing 
magnetization, could be detected while the temperature was still near its maximum 
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value: When the temperature was raised for a second time, only a few 
discontinuities occurred, in the same direction as before, as long as the 
temperature did not approach the previous maximum more closely than a few 
degrees centigrade. Eventually when the temperature had been raised and 
lowered several times over a given range, no further effect was observed in that 
range. 

It has been known for some time (Bush and Tebble, unpublished work) 
that if in the course of a slow traversal of the hysteresis loop the magnetizing 
field is allowed to fall by an amount whose magnitude varies with different 
materials, no Barkhausen pulses can be detected. If the specimen investigated 
was brought into a ‘reversible’ region of this kind, by lowering the applied field, 
and the temperature was then raised, no discontinuities were detected for a rise 
of several degrees centigrade. A further analogy between the effect of small 
changes of temperature and small changes of field was shown by the fact that, 
if the specimen was first cycled and then brought to the coercive region at a 
raised temperature such as 40°c, a decrease of temperature of as much as 20°C 
gave rise to hardly any detectable discontinuities at all. Any discontinuities 
observed were again in a direction of increasing magnetization. 

The second specimen investigated was of hard-drawn Swedish charcoal bnght 
iron, with a somewhat lower coercivity than the nickel. Similar results were 
obtained, except that a rise of about 100°c at remanence produced a few small 
discontinuities so distributed in direction as to decrease the magnetization of 
the specimen. This effect was not repeated when the temperature was raised a 
second time. 


§3. THEORETICAL DISCUSSION 


These results can be explained on the basis of two rather different mechanisms, 
one involving the thermal fluctuation field discussed by Néel (1950) in a 
treatment of the after effect in the Rayleigh region of ferromagnetics, and the 
other (suggested to the writer by Dr. P. Rhodes) arising from the decrease of the 
obstacles to domain wall movement, which usually accompanies a rise in 
temperature. 

The action of both these effects can be understood by considering the motion 
of a point representing the position of a domain boundary on an energy curve 
of the type described by Stoner (1950, p. 134). In a specimen subjected to a 
gradually increasing field, all the domain walls will occupy positions below points 
of inflection of the curve. As the position of each boundary is determined by an 
expression of the form H=const. dE/dx, where E is the energy of the system 
and «x the position of the domain boundary, and as Barkhausen jumps take place 
only when a boundary passes a point of inflection, it is clear that with a small 
decrease of field the effects should be reversible. This general picture is 
consistent with the experimental results, though it should be stressed that these 
give evidence of the absence of discontinuities above a certain minimum size, over 
a field range of up to 3 or 4 oersteds, and not of complete reversibility over that 
range. 

If the temperature of a specimen in a magnetic field is raised the thermal 
fluctuation field, which causes an oscillation of the domain boundaries, will be 
increased; in addition the slopes of the peaks of the energy curves for each boundary 
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will in general be decreased. For both of these reasons previously 
stable boundaries may execute irreversible jumps, though not necessarily 
instantaneously, in a direction of increasing magnetization. Neither of these 
mechanisms will of course operate for a fall of temperature, whether or not this 
was preceded by a rise in temperature at that particular field. 

If the applied field is lowered slightly before the temperature is raised, it is 
clear that a temperature rise large enough to cause the formerly most critically 
placed domain boundaries to jump, should now no longer be sufficient, since 
through the action of both the mechanisms mentioned above, all domain walls 
will have been moved to a more ‘ reversible’ region by the decrease of field. This 
again is in accordance with the experimental findings. 

At remanence, all domain walls will be in the troughs of their energy curves, 
and a much larger rise in the fluctuation field, or decrease in slope of the energy 
curves will be required to cause a Barkhausen jump. ‘The existence of a 
demagnetizing field will cause any jumps that do occur to be in a direction of 
demagnetization. 

A rough estimate of the relative importance of the two mechanisms discussed 
can be obtained by comparing the decrease of the coercivity with the increase 
of the internal fluctuation field with increase of temperature. Differentiation of 
the expression given by Néel (1950) for the fluctuation field leads to 
dH,/dT =(7k/vT)' where T is the absolute temperature, R is Boltzmann’s 
constant, and v will be taken as the volume of those Barkhausen discontinuities 
which make the largest contribution to the susceptibility. For both materials 
investigated wv will certainly be less than 10-!°cm? (Bush and 'Tebble 1948), 
which at room temperatures gives dH,/dT~10‘ oersteds/° c. 

Corner (1949) found that for nickel which had been placed under a strain of 
5 kg/mm? for several hours and then released, the rate of change of the coercivity 
with temperature is given by dH,,/dT~2 x 10-3 oersteds/° c at room temperature, 
and that for hard drawn iron dH,/dT~2 x 10-*oersteds/°c also at room 
temperature. 

These figures indicate that at least for the second material the boundary 
oscillation effect is probably of subsidiary importance. For other materials, 
however, such as permanent magnet materials with their necessarily smaller 
domain volumes, the thermal fluctuation effect may be of paramount importance. 


§4. RELEVANCE OF THE EXPERIMENTAL RESULTS TO OTHER WORK 


The experimental results described above, as well as the theoretical 
discussion, indicate that the changes of magnetization resulting from a decrease 
of temperature at any part of the hysteresis loop are likely to be almost entirely 
reversible, and that the irreversible component of the magnetization change will 
be a far smaller proportion of the total magnetization change than for a rise of 
temperature. This implies that the thermodynamic relations derived by Stoner 
and Rhodes (1949), between the magneto-caloric effect and the rate of change of 
magnetization with temperature at constant field, should be applicable over the 
whole region of the hysteresis loop instead of merely the saturation region, 
provided that the measurements of the magneto-caloric effect are carried out for 
small decrements of field, and provided that the measurements of (d//dT),, are 
carried out for small decrements of temperature. 
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ABSTRACT. The first set-up of a very simple high resolution 6-spectrograph is described- 
The instrument employs the focusing action of a magnetic prism with circular plane poles 
of 20 cm diameter. For this area of pole surface the dispersion is large, giving advantage to 
the width of source that may be used; and the instrument can focus electrons with energies 
up to4 Mev. It has been used with resolving powers Hp/A(Hp) up to 1300 as measured 
from microphotometer traces on the plates. It appears to be the only instrument that has 
been operated at these resolutions as a spectrograph. 

As a f-spectrometer, the resolving power should be 15% higher, but this has not yet 
been tested. As far as can be seen it has a luminosity roughly comparable with that of the 
ketron, but unlike that instrument it has considerable flexibility of experimental arrangement 
promising improved performance. 

It has been used to look for the fine structure of the K-conversion line corresponding to: 
the 1414 kev y-transition in the disintegration RaC.C’ reported by Latyshev and others. 
In agreement with the work of Bashilov and others who employed the ketron, and of 
Serebriskaya, Shpinel and Forafontov who used a spectrometer of completely novel design, 
this line was found to be highly monokinetic. 


§1. THE ORIGIN AND RELATIONS OF THE METHOD 
N the f-spectroscope described here a single magnetic prism with plane 
if circular pole faces is used to focus monokinetic f-rays from a line source 
placed in a direction parallel to the axis of the poles. 
Consideration of the field of two ideal circular poles separated by a distance 
equal to their radius (Siday 1947) showed that under appropriate conditions, 
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and as far as rays in the median plane are concerned, the quadratic aberration 
in the image could be eliminated and the cubic aberration reduced by a suitable 
choice of image position corresponding to unit magnification. The appropriate 
conditions were: 

(a) that the principal ray of the bundle focused should be initially directed 
towards the axis of symmetry of the field; (6) that a particular deviation of the 
principal ray of about 1-33 radians should be chosen. This means that with a 
given strength of magnetic field, the best focusing is obtained for a particular 
value of momentum, or Hp. 

These considerations formulate the essential requirements for the 
arrangement to function as a f-spectrometer of high resolution. A Geiger 
counter slit, or other line detector, is placed at the image position, and the: 
instrument is operated as a spectrometer by varying the magnetic field. 

Subsequent theoretical considerations (Jennings 1952) and experimental 
work with an electronic model apparatus (Ehrenberg and Jennings 1952) have 
demonstrated that ‘cubic’ focusing occurs under more general conditions. 
A new cubic focus can always be found by changing both the direction of the 
principal ray and its deviation. ‘Thus a series of images corrected for quadratic 
aberration exists, and if a photographic plate is placed to receive these the 
instrument so constructed will behave as a f-spectrograph; with a fixed field 
the spectrum of f-particles over a significant range of Hp is dispersed along 
the plate with approximately constant resolution. 

The elimination or reduction of an aberration in f-image formation appears 
to offer the f-spectroscopist the advantages of greater gathering power and 
greater resolution. Unfortunately, in spite of many important recent 
improvements in f-spectroscopy, modern instruments appear to fall into two 
categories offering one or the other advantage, but not both. ‘The aim underlying 
the development of the instrument described here has been to obtain high 
resolution; the gathering power is limited by the poor focusing properties 
(z-focusing) of the field on the trajectories not lying in the median plane. As an 
introduction and summary the present instrument is compared first with other 
instruments serving a similar purpose. 


(i) The Semi-Circular Focusing Spectrograph 

This instrument of classic importance is attractive because of its great 
simplicity of design, and because it gives a constant resolution over a very wide 
range of Hp. The essential disadvantage from which it suffers is the quadratic 
aberration of its image formation, making the line shapes unsuitable for the 
highest resolution. 

In its present form the new instrument gives cubic focusing, but the range 
of constant resolution is reduced to a 5°, variation of Hp, although the degradation 
of the image for greater variations is slow. Also it is worse than the semi-circular 
focusing instrument with respect to its 2-focusing properties. 

On the credit side, the present arrangement of two cylindrical poles of 
20cm diameter, in which no special trouble has been taken to avoid the effect 
of the magnetizing coils and yoke on the axial symmetry of the field, is even 
simpler than that required for the semi-circular focusing spectrograph where 
attention has to be paid to uniformity of field. 

Another attractive feature of the new instrument is the large separation of 
source and image which enables a modest magnet to focus high energy electrons 
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(4 ev in the present case). This, of course, is one of the factors giving rise to 
the high dispersion of the instrument. If As is the change of image position 
along the photographic plate measured in centimetres for a change A(Hp) in Hp, 
it is convenient to define the dispersion as | As/A(Hp)]Hp cm. For an instrument 
of given type the dispersion is directly proportional to the size (linear dimension) 
of the instrument. The dispersion of the semi-circular focusing spectrograph 
is roughly equal to the distance between source and image. ‘The new instrument 
with a 20cm pole diameter has a dispersion of approximately 33cm, so is 
equivalent in this respect to a very large semi-circular focusing spectrograph. 


(ii) High Resolution Spectrometers 


In the very considerable variety of recent development in instrumental design 
for B-spectroscopy the instrument described here appears to be unique in its 
function as a f-spectrograph. In its more limited use as a B-spectrometer it 
can be compared with other high resolution instruments. Resolution is 
conventionally defined as the fraction A(Hp)/Hp where A(Hp) is the width of a 
monokinetic spectrum line at half its height above background, and resolving 
power as the reciprocal of this quantity. Recent writers have been using the 
term ‘high resolution’ or ‘very good resolution’ to cover widely different values, 
so for definiteness it is taken here that an instrument does not fall into this class 
unless it has actually been used to give a resolution (in the presence of all 
attendant effects which cooperate in spreading the image) that is better than 0-2% 

The principal instrument at present achieving this resolution is that developed 
by Korsunsky, Kelman and Petrov (1944) and later by Dzhelepov and Bashilov 
(1950), and is called by them a ketron. It is a modified semi-circular focusing 
spectrometer in which the field in the median plane is given a (one-dimensional) 
gradient in a direction perpendicular to the line joining source and image. A total 
angle of collection of 30° in this plane is used, and theoretically an aberration-free 
image should be obtained. The shaping of the pole surfaces to give the 
appropriate theoretical field distribution, a feature common to a great deal of 
modern f-spectroscopic design, of course introduces a considerable complexity 
into the experimental set-up. 

It is interesting to note that, as in the instrument described here, the 
improved focusing in the median plane is achieved in the ketron only at the 
expense of a considerable deterioration of the z-focusing compared with that 
of the standard semi-circular focusing spectrometer. For this reason, both the 
ketron and the spectrometer described here are operated for high resolution 
with very small aperture angles in a plane perpendicular to the median plane 
(z-aperture). 

With principal aperture angles of 2° and 30° respectively the ketron gives a 
resolving power of 660. It appears that the practice with this instrument is to 
change the resolving power by changing only the z-aperture. 

With the new instrument using principal apertures of 1:5° and 6° respectively 
a best resolving power of 1300 was obtained as measured by microphotometer 
traces from the photographic plates. Taking count of the intensity—density 
characteristic of the emulsion used, this should correspond to a resolving power 
of 1500 in counter experiments. This last prediction has yet to be tested by 
detailed experiment. As regards dispersion, the ketron is essentially a 
semi-circular focusing instrument and will therefore have the small dispersion 
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associated with the latter. With poles 34cm x24cm, the distance between 
source and image used by Dzhelepov and Bashilov was 16cm, which we take 
as the dispersion of their instrument. For equal pole area the new instrument 
would have a dispersion of (34 x 24/m x 100)! x 33 =53, more than three times 
as large. On the other hand, at the same resolving power the ketron appears 
to have a higher luminosity than our instrument.* 

Quite recently some other experiments have been described in which a 
resolution of about 0-13°% was obtained. In achieving this, Shpinel (1946, 
1950) has developed a spectrometer of a novel design, but it is not yet possible 
to compare its performance with that of the ketron. 


(ui) The Logistics of the Experimental Arrangement Adopted 

It should be emphasized that the experimental arrangement described here 
must be regarded as a first set-up. It is explained below how the effective part 
of the field used can be regarded as agreeing with that given by two ideal plane 
poles separated by a distance equal to their radius. The choice of this field 
distribution depended on convenience rather than merit since, as was mentioned 
above, the behaviour of trajectories in it had been computed, and the knowledge 
so obtained facilitated the assembly of the instrument. A few preliminary 
experiments with very narrow angle bundles. sufficed to show that the instrument 
behaved very closely in accordance with the calculations. 

The calculations were made for bundles of $-rays with their principal rays 
directed initially towards the axis of symmetry of the field. ‘The distances of 
closest real approach to the axis of rays inclined at different angles to the principal 
ray were computed, and by placing stops in positions thus calculated angular 
apertures could be defined accurately. So it was again for reasons of convenience 
rather than merit that this direction for the principal rays was chosen in practice. 

It is not intended to discuss here how the performance of the instrument can 
be improved, except to point out that the cubic coefficient of aberration can be 
reduced and indeed eliminated for one particular Hp by bringing the poles closer 
together (Siday 1941, Jennings 1952), and that, by changing the direction of the 
principal ray in the sense that it would pass through more of the field, the 
z-focusing properties can be improved until for a single direction of a ray double 
focusing conditions can be attained (Ehrenberg and Jennings 1952)—although it 
may not be advisable to change the direction as much as this. ‘There does appear 
to be considerable flexibility in arrangement in the new instrument and we shall 
merely say that improvement in some directions is to be expected. 


§2. APPARATUS AND THE FIELD CALIBRATION 
(i) Description of the Magnet 

Figures 1 and 1 (a), show the magnet that was used as the magnetic 
prism. The pole pieces were of 10cm radius, of 10cm depth, and their 
separation was also 10cm. The poles were connected via the rectangular portions 
A and B to the yoke Y of diameter 12cm. 

The magnet was mounted so that its median principal plane was at right 
angles to the direction of the earth’s (total) field. 

The magnetizing coils consisted of three coils of eight, six, and four layers, 
the three coils fitting coaxially one over the other, and the innermost, of eight 


* These rather elaborate details of comparison are given as the Russian publications may not be 
readily available everywhere. 
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layers, fitting closely over the yoke of the magnet. The diameters of these coils 
were such that the resistance was about the same for each, of the order of 15 ohms. 
Each coil was wound over a tank through which cooling water was circulated. 


(ii) Magnetic Measurements and Adjustment of Field 

Standardization of search coils. An apparatus was devised whereby field 
measurements could be made in the median principal plane at 1 cm intervals 
along radial lines from the cylindrical axis of the magnet with sufficient accuracy 
for comparison with the computed field. 

A small search coil was used and its area found electrically by comparing the 
mutual inductance between the coil and a pair of Helmholtz coils, with that of a 
larger standard coil whose area had been directly determined from its dimensions 


Fig. 1. Design of prism magnet. Fig. 1(a). The spectrograph assembly. 


and number of turns. The Helmholtz coils and the standard coil were those 
used by Nettleton and Sugden (1939), and the method of calibrating the standard 
coil is described in detail by them. 

The magnetic field was measured by balancing the flux due to the removal 
of the search coil against the flux in the secondary of a standard mutual inductance, 
when a known current in the primary was broken. Using highly stretched 
‘elastic’, the search coil was catapulted from a fixed position in the field and at 
the same time the current was broken in the primary of the standard mutual 
inductance, the two charges resulting from these actions reaching their final 
values in less than j;second. Due to the inertia of the suspended system, the 
fluxmeter did not detect the differences in the rates of growth of the two charges, 
and provided the two charges were equal and opposite the fluxmeter recorded 
zero. No difference in result was obtained for large variations in the elastic 
tension employed, but a lower limit of tension was found, below which the two 
charges were seen to record independently. The tension employed in the elastic 
was always well above this limit. 
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The field in oersteds is given by the relation H=100Mi/A, where A cm? is 
the effective area of the search coil, M is the mutual inductance in wu, and 
2 is the current in amperes in the primary of the mutual inductance. An 
accuracy of field measurement of 1 part of 500 was satisfactory for our purposes 
and was achieved without particular difficulty. 

As the search coil had to be removed from the field with considerable 
acceleration, its path had to be guided in some way. ‘This was done by sinking 
the coil in a block of wood 3 inches long and 1 inch square in cross section, and 
shooting the block along a hollow brass tube of the same (internal) cross section. 
‘To one end of the block was attached the elastic and to the other end was screwed 
a brass cylinder with a central hole with 3in. Whitworth thread. Attached to 
the bottom of the brass guide tube and at right angles to it were three rectangular 
bars of $ inch cross section. These fitted into rails on a horizontal table, part of 
which extended between the poles of the magnet. The height of the table could 
be adjusted. With this arrangement the brass tube could be slid across the 
table in a direction parallel to the vertical* plane of symmetry of the magnet 
and at right angles to the length of the tube; the table height could be adjusted 
so that the central plane of the search coil lay in the median plane of the magnet, 


Fig. 2. Catapult for measuring the magnetic field. (Not to scale; the length of the tube is 6 feet.) 


The position of the tube relative to the magnet is shown in fig. 2 when the 
axis of the tube passes through the central vertical axis of the magnet. ‘lhe 
top of the tube was removed over the portion AB and the centre of the sides 
slotted over the same distance, while a slit 1 mm wide was milled in the centre 
of the bottom of the tube. A square piece of paper on which were ruled two 
lines at right angles—x =O in the direction of the plane of vertical symmetry, 
y=0 in the direction at right angles to this—was placed under the tube so that 
the line y=0 was visible through the 1 mm slit over its whole length. ‘The lines 
were marked off in 1cm intervals. K is a brass block similar to the block C 
containing the search coil, and it can be fixed rigidly in any position by means 
of the wing nuts W. A half-inch brass rod can be screwed into the block K at 
one end and the block C at its other end. ‘Thus the block C can be held rigidly 
in any position by clamping block K. 

The block K was moved until, on looking downwards through the j inch hole H 
drilled in the } inch brass rod S, the pair of cross wires above and below the hole 
were coincident with a desired mark on the paper. ‘The distance from the centre 
of the search coil to the centre of the crosswires was made equal to the distance 
from the centre of the magnetic poles to the origin «=0, y=0 on the paper, 


* The field plotting was done before the magnet was mounted with its axis in the direction of 


the earth’s field. 
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then if the crosswires were set on a point’y =0, « =7, it was known that the centre 
of the search coil was rcm from the centre of the pole. The brass rod was rotated 
by means of a special key, so that it unscrewed from the block C and screwed 
into block K. In this way the search coil was released. 

At the same instant a switch was automatically broken in the circuit of the 
primary of the mutual inductance and the secondary flux balanced as previously 
explained. Values of the field H were obtained for corresponding values of the 
radial distance 7 along the line y=0 up to r=32 cm. 

A small correction for the error in not placing the origin at exactly the correct 
distance from the centre of the magnet was found electrically by shooting the 
coil from a particular position x=7’, y=0, as read on the paper, and then from 
a position x= —r’, y=0. From the difference in the two readings of the magnetic 
field H, and a knowledge of the rate of change of H with 7 at this point 7’, the 
error x could be found from «=4AH/(dH/dr),.. 

A value r’ =9cm was chosen, for, in this region, (dH/dr),, is nearly constant. 
This correction was additive to all values of x as read on the paper and amounted 
to 0-4.cm. 

By moving the whole brass tube at right angles to itself and shooting the 
coil always from «=0 the variation of field in the y direction was explored. 
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As explained above it was decided to accept for preliminary work any field 
approximating reasonably well to the computed one. In fig. 3, curve I shows 
the field shape as computed, and curve II shows that determined experimentally 
along the y=0 axis. The field shape at right angles along the x =0 axis did not 
differ significantly from that of curve II, so that for the purpose of comparing the 
experimental and computed fields the asymmetry need not be further considered. 
It was noticed that, if the abscissae of curve I were multiplied by a constant 
factor 1-35, a fair fit with curve II was obtained (curve III). This means that 
the magnetic field shape in the plane of symmetry of the actual magnet is more 
nearly equivalent to the computed field shape for a magnet of pole radius 13-5 cm 
and separation 13-5cm than of 10cm. Therefore, the focusing properties of the 
actual magnet should be similar to those theoretically considered, provided all 
distances concerning electron trajectories are multiplied by a factor 1-35. 
It will be seen that the most significant disagreement between curves III and II 
occurs in the region between r=0 and r=10cm. However, the focusing 
properties of the prism have been considered only for electrons leaving the 
source at an angle in the plane of symmetry sin’! 0-3. The nearest approach of 
these electrons to the vertical axis of the magnet, according to computation, is 
4-27cm for the optimum Hp value, and 3-75cm for the greatest Hp value of 
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electrons to be received on the photographic plate. If these values of r are 
multiplied by 1-35 it will appear that the field of the actual magnet in the region 
r<5-4cm has no bearing on the performance of the instrument. This being the 
case (so long as the angular spread of the accepted sheaf of electrons is limited 
as described above) comparison of ‘scaled computed’ and experimental fields 
is best made by adjusting these to equality at r=5-4cm. This is done in fig. 4, 
where remaining discrepancies are very small indeed. 

Subsidiary experiments showed that the maximum change 1 in effective field 
for 1°c rise in yoke temperature was less than one part in 30000. Therefore, a 
hypothetical resolution as good as 1 part in 10000 should not be impaired by 
temperature effects if it is supposed that by water cooling the coils the temperature 
of the magnet can be kept constant to one degree centigrade for several hours. 

In order to maintain a reproducible relationship between magnetizing current 
and field, the magnet was first saturated by passing a total current of about 
45 amperes for a period of about half a minute through the three coils connected 
in parallel. This current, obtained from the 230 volt d.c. mains, was then rapidly 
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reduced to zero. It was found that a residual field of 14-7 oersteds was always 
obtained. Using the ‘ascending’ portion of the magnetization cycle it was then 
found that any required field could be set to at least 1 part in two or three hundred. 

The inside magnetizing coil acting alone was capable of producing a field at 
the centre of 850 oersteds, but for fields above this value all three magnetizing coils 
were used in parallel. By means of simple experiments it was found possible to 
attribute relative magnetizing powers K and L to the middle and outer coils 
respectively. A current 7, in the middle coil was thus equivalent to a current Kz, 
in the inner coil, and z, in the outer coil was equivalent to Li, in the inside coil. 
When all three coils were used the effective magnetizing current was taken as 
1, + Ki, + Lg. 

A graph showing the relationship between effective current and field is shown 
in fig. 5. ‘The magnetizing current elec to focus some of the well-known 
lines of Th(B + C) is indicated. 

The field was stabilized by keeping the current through the magnetizing coils 
constant during the time electrons were being recorded. ‘This current was 
supplied by a hundred-volt bank of accumulators. The current was measured 
by measuring the voltage drop across a standard resistance with a ‘Tinsley vernier 
potentiometer used in conjunction with a low-resistance high-sensitivity 
galvanometer. A change in current of 1 part in 40000 was readily detectable. 
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By connecting a high resistance in parallel with the control resistance and by 
using a gear reducing device for altering the former by small amounts, a 
magnetizing current of about half an ampere could be kept constant to 
considerably better than 1 part in 40000 for several hours. The magnetizing 
current could easily be altered in measurable steps amounting to 1 part in 20000 
of its average value. The maximum energy of electrons that could be focused, 
using a magnetizing current supply from one hundred volts, was 4-2 Mev. 
(iii) Spectrograph Chamber 

The vacuum chamber was constructed of $in. thick sheet brass, and parts 
were soft soldered. The external height of the chamber was 2:5in. The top and 
bottom plates of the chamber were of identical construction, as shown in fig. 6. 
Slots in. deep and 2in. wide were milled in these plates and the source and 
detector of electrons could be placed in any required position, relative to the 
axis of the magnet, by sliding them along the slots. ‘The chamber was lined 
with aluminium #in. thick, which could be removed by sliding it along the outside 
slots P in fig. 6. The lining of the sides extended along their whole length, and 


Fig. 7. Rays in vacuum chamber. 


the lining of the top and bottom plates was attached to these sidepieces. ‘I’o avoid 
covering up the slots Q and R (fig. 6) where they were required to hold the source 
and detector, the top and bottom plates were lined only in the region between the 
source and detector. 

Two coaxial brass rings S were attached to the outside of the chamber, one 
to the bottom plate and one to the top plate. ‘They were of 10cm radius and 
1-5in. depth, and the pole pieces of the magnet fitted closely into them. Thus the 
position of the chamber was fixed relative to the axis of the magnet. he vertical 
plane of symmetry AB (fig. 7) of the chamber was then coincident with the 
vertical plane of symmetry of the magnet. Removable brass packing pieces 
were placed between the chamber and the magnetic poles so that the median 
principal plane of the magnet coincided exactly with that of the vacuum chamber. 

Each end of the chamber terminated in a rectangular flange, and a rubber 
gasket of elliptical cross section fitted into a rectangular slot in each flange. An 
end plate could be clamped to each flange and the rubber gasket compressed to 
produce a vacuum seal. 

The pumping system consisted of a two-stage Metropolitan-Vickers type O3B 
oil diffusion pump joined to the vacuum chamber by a pipeline less than two inches 
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long. ‘The pump was of all brass construction and could be connected directly 
to the chamber without disturbing the distribution of the stray field of the magnet. 
The backing pump, however, was made of magnetic material and was placed 
at a distance of about six feet from the spectrometer. It was joined to the backing 
side of the diffusion pump by a long length of brass tubing. A pressure of less 
than 10-°mm could be maintained in the spectrograph chamber. 

In accordance with the plan of following the computations the angle of the 
chamber was made 100°. 

For the optimum Hp value, the computations correlated the values of the 
initial inclination of a ray to the principal ray with the nearest real approach 7,, 
of the electrons to the axis of the magnet. 7,, lies in the vertical plane of 
symmetry AB (fig. 7) of the chamber. Within the range studied, any required 
angle of emergence could be selected by placing an aperture in the plane AB 
and choosing its limits 7,,. The aperture was carried by an aluminium plate 
attached to the removable aluminium lining that fitted into the slots P (fig. 6) 
and was adjustable in the horizontal and vertical directions. Regarding the 
vertical direction, the computations show that most of the defocusing action of 
the prism takes place after the electrons have traversed the vertical plane of 
symmetry AB. Therefore the height of the aperture in the x direction was a 
good measure of the maximum angle of divergence of transmitted electrons, the 
angle being measured at right angles to the median principal plane. The edges 
of the aperture were bevelled to reduce scattering of electrons. 

N was an aluminium block that prevented secondary electrons, produced by 
y-rays striking the sides of the chamber, from reaching the detector P. Lead 
blocks L were placed between the source and detector. 

The source of electrons. ‘The performance of the instrument was investigated 
using the F line of thorium B. In a first series of experiments active deposit 
from radio-thorium was obtained on a platinum wire 2cm long and 0:1mm 
diameter. For more accurate work a second source holder was constructed so 
that a strip of platinum foil 0-1mm thick, 2cm long and 1cem deep could be 
activated along the 2cm edge. Only the edge of the foil was exposed to the 
radio-thorium, the contamination of the sides being prevented by compressing 
the foil between two thin pieces of rubber stuck to a brass plug. When the activity 
of the foil was very weak it was placed for a few seconds on a photographic 
emulsion. ‘The blackening of the emulsion had a sharp edge corresponding to 
the edge of the foil and this showed that the side of the foil was not greatly 
activated. 

After activation the foil was removed and mounted on a brass plug. It could 
always be replaced on this plug in an exact relative position by compressing it 
between two flat brass surfaces. In this way a non-contaminated source mounting 
could be placed in the spectrometer. ‘The plug could be placed on a holder which 
slid into the spectrometer vacuum chamber along grooves Q (fig. 6). Its position 
relative to the axis of the magnet could be adjusted by a micrometer screw that 
could be attached to the end flange of the chamber before screwing on the end 
plate. A slot milled in the base of the micrometer plate allowed the plate to 
come into flat contact with the brass flange F without compressing the rubber 
gasket. ‘The micrometer could be set to 0-01 mm. 

By means of small adjusting screws the source could be made accurately 
parallel with the vertical component of the magnetic field. 
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§3. PERFORMANCE OF THE INSTRUMENT USED AS A SPECTROGRAPH 
Photographic emulsions were used to record the focused electrons. Ilford 
X-ray plates and Ilford G-5 Nuclear Emulsion plates were used with the 
spectrometer. The former were very sensitive but had a large grain size. The 
nuclear emulsion plates were much slower, but had a very small grain size and 
an emulsion thickness of 50. For detection of weak lines against a continuous 
B-spectrum by eye, the X-ray plate was only slightly inferior, but the fluctuations 
due to grain size made photometric measurements of very weak lines difficult. 

The plates were 3}in. long and lin. high. They fitted into a plateholder 
in a camera which was loaded in a darkroom. ‘The camera was transferred to 
the vacuum chamber and its positions could be adjusted by a micrometer screw 
reading to 001mm. Once a suitable vacuum had been obtained the front of 
the camera could be rapidly opened, by rotating a rod passing through a Wilson 
vacuum seal, and the plate exposed to the electrons. After exposure the camera 
could be closed and removed to the dark room for unloading. For a pole radius 
of 10cm the focal length for best focusing for a principal ray was calculated to be 
10-8 cm and therefore, taking account of the experimental field distribution, the 
soutce was placed at a distance 210-8 x 1-35=29cm from the axis of the 
magnet. : 

The investigation of the field shape enabled us to place the photographic 
plate in approximately the correct position, at the same distance from the axis 
of the magnet as the source. The angle of emission of electrons in the median 
principal plane was adjusted to 6° and several separate exposures were taken for 
different positions of the camera. The plates were microphotometered and it 
was found that the optimum position of focus differed from the theoretical by a 
few millimetres. The angle of emission was now increased to 18° and the focusing 
procedure repeated. The optimum focus for the larger aperture coincided with 
that for the smaller, but the image was wider. 

The focusing field was increased by 10° and the camera refocused. The 
plate-holder was then placed so that the line joining the foci corresponding to. 
Hy and 1-1 Hy was coincident with the emulsion. With the plate at this angle, 
the same camera position gave the best focus for a magnetizing field 10% below 
the optimum. One photographic plate could thus record a spectrum range of 
20% in momentum. The angle of incidence of the central ray at the plate in this 
arrangement was about 46°. ‘The angle of deviation of the principal ray, 
characterized by b = 0, for optimum focusing was in agreement with that computed. 

The photographic plates were photometered with a Kipp automatic recording 
microphotometer. ‘Iwo magnifications were possible, 50 times the distribution 
on the plate and 5 times the distribution on the plate. Horizontal and vertical 
scale lines were put on the recording film by two lamps; one of them flashed 
intermittently, the flashes being controlled by the motor. On the larger 
magnification the separation of the vertical distance lines was 5 mm corresponding 
to 0-1 mm of plate traversed. On the smaller magnification the separation was 
2mm corresponding to 0-4 mm of plate traversed. 

‘Tests were made to ensure that the shape of the trace was not affected by small 
alterations in the width and length of the analysing slit of the microphotometer. 
Usually a slit width of about 10 and length 0:2mm was employed. As was 
to be expected, the focused line varied in shape and intensity from point to point 
along its length (the z-direction). Two marked discontinuities in intensity were 
found at distances from the centre of the line or the median principal plane, which 
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depended upon the source-length and the angular aperture in the z-direction. 
The lines were not straight but slightly curved. 

Reduced copies of microphotometer traces of the F line are shown in 
figs. 8 (a), (6) and (c). The source length was 2cm and angle of emergence in 
the median principal plane was 6°. The total angle of emission in the z-direction 
was 1-5°. (a) is the trace across the centre of the line, (b) is 2mm from the centre. 
(>) has the same intensity above background as (a), but the line shape of (b) is 
sharper, particularly on the low energy side. (c) is the trace measured at 4mm 


from the centre. 
fol (c) 


(a) i (b) 


Fig. 8. Reduced microphotometer traces of F line. 


b 
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Fig. 9. Reduced traces of J,, J, and Ja lines. 


The width of the trace lines at half their height above background will be 
called the trace-breadth. ‘This will be different from the true half width of the 
line and will vary with the exposure, even over the linear portion of the density— 
log (exposure) characteristic. 

The dispersion factor was obtained by altering the magnetic field and 
measuring the line shift, and has also been obtained by recording two B-ray lines 
whose absolute energies are known and measuring their separation on the trace. 
The agreement between the two determinations was good. ‘The J and Ja lines 
of thorium B had the most suitable energy difference, and the trace obtained 
is shown in fig. 9. The difference in energy (Ellis 1932) is 3kev and the 
separation of their peaks is 2.5mm on the photographic plate, giving a value 
D=330mm. Therefore each division on the trace (corresponding to 0:-1mm 
on the photographic plate) corresponds to Hp/A(Hp) =3 330. Figure 8 (6) has a 
trace breadth corresponding to a resolving power, measured as Hp/A(Hp), of 
approximately 1 300. 

A large number of exposures were taken for various source lengths and 
aperture sizes. The general results are summarized below. 

Trace-breadths corresponding to better resolving power than about 1300 
could not be obtained by making the angle of emission in the median principal 
plane less than 6°. 
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The variation of angle of emission in the z direction did not change the 
intensity of the line in the median principal plane by a large amount until the 
angle of emission was greater than 2°, after which the resolving power decreased. 

The decrease in resolving power with angle of emergence in the median 
principal plane was not very rapid up to an angle of 12°. Thus with an angle of 
emergence of 12° a resolving power of 800 was obtained. For an angle of 
emergence of 18° the resolving power, measured from the half-width, was 350, 
but the line had a long tail on both the high and low energy sides of the peak. 

For changes in field of +3°% about the optimum Hp, there was no change 
in resolving power. For larger changes the resolving power decreased and the 
lines became asymmetric, having a tail on the high energy side. This tail could 
be considerably reduced by asymmetrical adjustment of the aperture in AB, 
that is by choosing a suitable direction of the principal ray, but this disturbed 
the focusing for fields of the order of H). 

Figure 10 shows the trace obtained after a double exposure to the F line. 
The field was changed by 1 part in 1 200 and the lines are clearly resolved. 
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Fig. 10. Reduced microphotometer trace; Fig. 11. Geiger counter determination 
resolution test on the F line. of the thorium F line. 


If a radioactive material has a short half-life it may be necessary to use several 
sources with the same photographic plate. It is then necessary to be able to 
replace the source exactly. ‘T'o test the ability to do this, an exposure was taken 
using the Thorium F line. ‘The camera was closed, the pumps were switched off 
and air admittted to the chamber. ‘The source was completely dismantled and 
reassembled and a second exposure of the same length as the first one was taken. 
The microphotometer trace showed that the source was replaceable with such 
accuracy that no change in line shape could be detected. 


§4. GEIGER COUNTER MEASUREMENTS 
Geiger counters were also used for recording the focused electrons. Two 


counters were connected in coincidence so that only particles that passed through 
both counters were recorded. 


‘The plate-holder in the camera could be removed and the two counters fixed 
in the camera. ‘The counters were cylindrical in shape, the axis of each counter 
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being parallel to the axial component of the magnetic field. The counters were 
made of copper and were 1-25in. long and in. in diameter, with Ain. wall 
thickness. ‘The central wire was insulated from the walls by a small metal—glass 
seal and a small glass bead was attached to its other end. 

Mica windows were fitted to the counters. The thinnest window that could 
be fitted to the first counter without breaking had an average thickness of 
1-4 mg/cm? and the thinnest fitted to the second counter an average thickness of 
2:9mg/cm?. ‘The counters were pumped out for long periods and filled with 
argon and about 10° of alcohol vapour. ‘They were filled to a pressure of 
about 10cm. 

Two separate power units supplied the potentials for the central wires. They 
were standard A.E.R.E. power supply units type 1082 A. The scaler employed 
was also a standard A.E.R.E. scale of 100. 

The scanning slit in front of the first counter was 40 1 wide in brass. It was 
placed at the focus of the electrons and was bevelled on the side furthest from 
the source. The permitted angle of acceptance of electrons by the slit was 
limited to about 6°. This was considered necessary because high energy 
electrons have a considerable range in brass. The slit was about 2cm long and 
parabolic in shape. It was very difficult to be sure that the curvature of the 
slit was exactly the same as that of the natural line. The line curvature will be 
different from the photographed curvature for in the latter case the electrons 
are incident at an angle with the plate but strike the slit normally. The slit was 
adjustable in all directions. 

The Geiger counters were not used to make a comprehensive study of the 
properties of the prism, but a line obtained with counters is shown in fig. 11. 
The magnetizing current was changed in steps of 5 x 10-° amp near the peak 
and in steps of 2 x 10-4 amp at the base of the line. 

The line has a half-width corresponding to Hp/A(Hp)=900. A piece of 
brass 2mm wide was then placed across the centre of the slit, so that the rays 
focusing in the median principal plane were prevented from reaching the 
counters. ‘The counting rate was decreased by about 20°, but the resolving 
power was increased to 1100. Converting microphotometer traces to relative 
intensity units, using the measured characteristics of the blackening of the 
photographic plates, showed that under best conditions the true half-width as 
distinct from the trace-breadth corresponded to a resolving power of 1500. 
It seems likely that by more careful shaping of the scanning slit a resolving 
power nearer this figure could be obtained with counters. 

With counters, for the F line of thorium B, the ratio of the height of the peak 
to the height of the continuous spectrum was found to be about 300. 


§5. THE SEARCH FOR FINE STRUCTURE IN THE INTERNAL 
CONVERSION LINES OF RADIUM C 
During 1948-49 Latyshev and his co-workers published a number of papers 
claiming that many of the well-known internal conversion lines of radium C.C’ 
exhibited a fine structure, the component lines being equally spaced. Latyshev 
has now (1950) withdrawn this claim and suggests that the multiple sharp peaks 
found experimentally were due to insufhicient counting and to fluctuations in 
their magnetic field. Before this withdrawal preliminary investigations had 
already been made using the new prism spectrometer, and a brief report on our 
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definitive findings regarding the structure of the 1414kev and 606 kev internal 
conversion lines is given below. 

Because of its longer lifetime radon is more convenient to use as a source 
than is a deposit of radium C, though the straggling of the electrons in passing 
through the walls of the container of the radon gas will produce a broadening of 
the image. A thin-walled glass ampoule was therefore filled with 100 millicuries 
of radon. The wall thickness was 7 x 10-4cm, or about 3mg/cm?. ‘The tube 
was made by drawing out thin-walled soda glass tubing. Many tubes were drawn 
and examined under a microscope. A rough check of the wall thickness of the 
tube was obtained by measuring the maximum range in air of the «-particles 
emerging from the tube. The length of the ampoule was about 2cm and its 
diameter 0-012cm. The fine glass tube was mounted on the same source holder 
and in the same position as the platinum wire that was used in the focusing 
experiments. It was held in place by tap grease. 

The K, L and M conversion lines of the 1414 kev y-transition were recorded 
using Ilford G5 Nuclear plates. Several exposures were taken with this source 
and a microphotometer trace is shown in fig. 12. The trace breadth of the 
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Fig. 12. Reduced microphotometer trace of the conversion lines of the 1 414 kev y-transition RaC.C’ 


K conversion line is about 10kev. If, however, the trace 1s converted into 
relative intensity units, the line shape is considerably altered and the half-width 
is only 5kev. From the results of White and Millington (1928), a line of zero 
width should be broadened to a half-width of about 4kev for electrons of this 
energy passing through a ‘light’ material of thickness 3mg/cm?. The 
instrumental half-width corresponding to a resolving power of 1000 should be 
about 1:8kev. Therefore the line had about the half-width to be expected from 
a single line for electrons passing through the walls of the glass ampoule. 

To avoid the line broadening due to the straggling of the electrons in the 
giass, it was decided to activate the edge of a piece of platinum foil. The edge was 
exposed to radon gas and a negative potential of 1000 volts was applied to the 
foil; 200 millicuries of radon were used and the period of activation was one hour. 
The y-ray activity of the active deposit on the foil was equivalent to that of about 
25 mg radium in equilibrium with its products when first removed from the radon. 

The Geiger counters were used in coincidence and the timing of the readings 
carefully noted in order to correct for source decay during counting. 

With the K conversion line of the 1414 kev y-transition the number of counts 
taken was such that the statistical mean error at the background level was 
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approximately 6°%, and at the peak of the line 3%. Values between these were 
obtained at other points. The ratio of the height of the line to the height of the 
continuous spectrum was 6:9:1. ‘This ratio is a measure of the resolving power 
of an instrument. ‘The value obtained by Latyshev was 2:1. Figure 13 shows 
the line obtained (curve II); curve I represents Latyshev’s line. The continuous 
spectrum has been deducted and the curves normalized to have the same peak 
intensity. The half-width of our line corresponds to a resolving power of 
approximately 800. 
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Figs. 13 and 14. Geiger counter determinations of the K conversion lines for the 1414 kev 
and 606 kev y-transitions of RaC.C’. In each case the outer curve is that due to Latyshev. 


The platinum foil was reactivated and the K conversion electrons from the 
606 key y-ray were examined. Figure 14 shows the result in curve II, after 
subtraction of the continuous spectrum and normalizing the line, so that its 
peak coincides with the line obtained by Latyshev. ‘The number of counts 
taken gave a statistical mean error at the background level of approximately 5% 
and at the peak 3%, and values between these at other points. ‘lhe half-width 
of the line corresponds to a resolving power of approximately 1000. 

Bashilov, Antoneva and Zolotavin (1950) state that they have examined the 
606, 1120, 1414 and 1760 kev K conversion lines of radium C at a resolving 
power of about 400. ‘They found no fine structure. 

Serebriskaya, Forafontov and Shpinel (1950) also investigated the same 
lines but because they used radon gas in a glass ampoule their resolving power 
was inferior to that of Bashilov, Antoneva and Zolotavin. 
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The Re-investigation of some Photoelectron Lines 
of Thorium (B+C) 
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1st Fanuary 1952 


ABSTRACT. The $-ray spectrum of thorium (B+ C) has been re-investigated using the 
f-ray spectrograph described by Siday and Silverston. The range 200 kev—300 kev was 
studied in detail using photographic recording and a previously unobserved y-ray of energy 
294-6 kev has been attributed to the transition thorium C.C”. The conversion lines of the 
yF¥ ray of thorium B.C in the Ly, Lyy electronic shells have been resolved and their relative 
intensities have been found using a Geiger counter. The conversion of the same y-ray in 
the M,, My electron shells has been measured and the lines have been partially resolved. 
The M photoelectron line associated with the transition thorium C’.D has been examined, 
and its suspected doublet structure confirmed. 


HE f-spectrograph described by Siday and Silverston (1952) has been 
| used to record the lines in the region 200 to 300 kev. The thorium-active 
deposit was carried on a fine platinum wire of diameter 0:l!mm. The 
relative energies of the lines are given in the table and are compared with the 
values given by Ellis (1932 a) and Surugue (1937). Some microphotometer 
traces are shown in fig. 1. Not all the lines that are indicated could be identified 
from these traces alone, but a large number of traces at ten times this 
magnification were studied and the ordinates were averaged. The nomenclature 
is similar to that devised by Ellis (1932 b) and, except for the addition of a new 
line, is exactly the same as that adopted by Surugue (1937). 

The difference in energy between the I and J lines, as given by Ellis (1932 a), 
was used to measure the dispersion of the spectrograph, and the energies of all 
the other lines were obtained on the assumption that the I and J lines differed 
in energy by 12-30 kev. This is the difference in energy between the L; and M, 
electron shell levels for element 83. The absolute energy of the I line was 
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Photoelectron Lines of Thorium (B+C) 


Energy (kev) 
Name Silverston Ellis Surugue 

202-1 201-7 201-9 
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Energy (kev) 
Name Silverston Ellis 


Jas 241-4 
Jag 255-4 
Ja, 257°7 
Jb 260-9 
Jb, (E) 263-9 
Jb, 271-7 
Jb, 272-6 
Jb, 275°7 
Jb, 279-9 
Je 283-7 
Je, 284-3 
i 295°5 
i 300-3 

221-7 kev 
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241°5 
25920 
PC 
260-9 
263-6 


2724 
282°8 


YET] 


Surugue 


241°9 
256-0 
258°4 
261-0 
Pipes, 
Ped 
275°8 
278:°4 
283°8 
284-6 
296-1 
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Fig. 1. 


Microphotometer traces of part of the B-spectrum of thorium (B-+C). 
be found in the position marked H, and the presence of a line in the position marked J* 


was not convincing.} 


x3 


No line could 


+ Note: In the lower two tracings the J on the right-hand side refers to J*. 
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assumed to be 221-7kev, this being the figure given by Ellis (1932 a). Except — 
in one or two instances the relative intensities of the lines have not been measured. ~ 
The chief results of the new survey are as follows: 

Ellis found the line Ia, to be of energy 233-2 kev and Surugue (1937) gave a 
value 229-0kev. The energy discrepancy between these two lines is so large 
that, although both authors have given them the same name, Ia,, it would appear 
likely that this refers to two different lines. Neither of these was found, however, 
but a line, which has also been referred to as Ia,, was found at 231-0kev. Its 
origin has not been explained. 

The line listed by Ellis (1932 a) at 267-2 kev could not be found, nor could 
the H, line listed by Ellis and Surugue at about 213 kev. » 

Surugue interpreted the Ja and Ja, lines as being the N and O electron shell 
conversion lines of the yF ray of element 83. He gave the ratio of their intensities 
as greater than 18:1. The ratio has been found to be less than 4:1. 


The Ha Line 

There appears to be a weak line on the high energy side of the H line. Its 
energy is estimated to be 0-6 kev higher, that is 209-5kev. If this line were due 
to the conversion in the K shell of the y-ray from element 81 the energy of the 
y-ray would be 294-6kev. The L; conversion of this y-ray would give a line 
at 279:3kev. The energy of the line Jb, was found to be 279-9kev, in good 
agreement with this predicted value. The origin of the Jb, line has not been 
previously explained. 

Measurement of the energies of the «-particles emitted in the transformation 
thorium C.C” by several authors suggests that there should be a difference in 
the nuclear levels corresponding to 294kev. Ellis (1932 a) suggested that the 
ray yH, of 298kev could be associated with this energy difference of the 
a-particle groups. ‘This is not very convincing for the following reason. He 
tabulated eight different y-ray energies and compared them with the energy 
differences deduced from «-particle measurements. All the y-rays with the 
exception of yH, had an energy a few kev less than the corresponding energies 
deduced from «-particle measurements. In the case of H, the energy was greater 
by about 4kev. Furthermore, Surugue has attributed the H, f-ray to a 
thorium C”.D transition. 

The H photoelectron line is due to the K conversion of the yH ray of 
element 83 and its L; conversion electrons account for most of the intensity of 
the Jc line. The ratio of the intensity of H: Jc is given as 10:1 by Surugue and 
about 8:1 by Ellis. ‘These ratios are rather greater than the usual ratio of 
intensities for K: L conversion at this energy, but the fact that the Ha line was 
not previously resolved from H may account for the discrepancy. 


The J Line 

The microphotometer trace of this line is shown in fig. 2, and fig. 3 shows 
the form of this line after conversion to relative intensity units. The broken 
curves are the instrumental line shapes, and their intensities and separation 
have been adjusted so that the sum of their ordinates gives the experimental 
(full line) curve. ‘The J line is due to the M; and M,, conversion electrons of 
the yF ray of thorium C. The separation of the dotted curves in fig. 3 is 
0-25 kev. From x-ray data the separation of the M; and M,, levels is 0-3 kev. 
The ratio of the intensities of the two components is 4:3: 1. 
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Fig, 2. Microphotometer traces of the J and Ja lines. x 30 
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Fig. 3. The J line after conversion to relative intensity units. 
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Fig. 4. Microphotometer trace of the I and Ia lines showing the Li, Lr levels resolved. 30 
Note: Ly, refers to the Ln level. 
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The I Line 


The I line is due to the L; conversion of the yF ray from Thorium C. The 
Ly conversion level is at 0-6kev higher energy, the separation of the levels 
requiring a resolving power of about 600 expressed as Hp/A(Hp). 

The microphotometer trace is shown in fig. 4. 

The lines were also measured with two Geiger counters in coincidence. 
The curve obtained after deduction of the continuous f-spectrum is shown in 
fig. 5. The number of counts taken on the high energy side gave a statistical 
mean error of less than 2° at the continuous background level. 

The analysis of the line into the two components (broken curves), gives 
exactly the correct separation of the levels. The small dotted curve has been 
drawn to have the same shape as the larger curve but its area is one-eighteenth of 

the area of the larger curve. 
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Fig. 5. Photoelectron spectrum of the I and Ia Fig.6. The M and L lines. The weaker 
lines after deduction of the continuous (dotted) curve on the low energy side 
B-spectrum. The statistical mean error is has been drawn the same shape as the 
less than 2% at the background level. The Lline, but with one-fifth of its inten- 
analysis of the line into the two components sity. The addition of the ordinates 
(dotted curves) gives exactly the separation approximates closely to the M line. 


of the Li, Lu levels. Their ratio of inten- 
sities is 18: 1. 


The ratio of the intensities of the L; and Lj, conversion lines was thus found 
to be 18:1. Forafontov and Shpinel (1950) have found this ratio to be 
approximately 20:1. ‘The values obtained photographically by Ellis (1932 a), 
Surugue (1937) and Arnoult (1939) were 15-6, 16 and 11 respectively. 


The M Line 


In order to explain the fact that in the transition thorium C”.D the L, 
M and N photoelectron lines are equally spaced and appear to differ in energy 
by an amount equal to the difference in energy between the K and L, electron 
shells of thorium D, it has been assumed that two y-rays (yL and yM) differ in 
energy (fortuitously) by this amount. On this assumption the M line is of double 
origin, and its intensity relative to the L line (roughly equal according to Ellis 
and Surugue, 1: 1-3 according to Flammersfeld (1939) and 1: 0-83 according to. 
Martin and Richardson (1948) is consonant with this view. Only about 
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one sixth of this intensity is likely to be due to the L; photoelectrons from yL, 
the remainder arising from the K photoelectrons from yM. 

In order to check this assumption directly, an attempt was made to resolve 
the M line photographically. Three exposures were taken and microphotometered 
at several points along each line. Three exposures were also taken of the L line. 
The traces have been averaged and converted to relative intensity units for each 
line, although the relative intensities of the L and M lines have not been compared. 

Figure 6 shows the two curves. The weaker (dotted) line on the low-energy 
side has been drawn the same shape as the L line, but with one fifth of its intensity. 
If the, ordinates of this line are added to those of the L line a line very similar 
to the experimental M line is obtained. 

The bump at the base of the line on the high-energy ede was not interpreted. 

If fig. 6 represents the correct analysis, the separation of the L,; conversion 
line of the yL ray and the K conversion line of the yM ray is approximately 
0-43 kev. 

If the energy of the M photoelectron line is 494-17 kev (Ellis 1932 a), then 
the L; conversion line of the yL ray will be 493-74kev. ‘Taking the K and L, 
absorption energies as 87:50 kev and 15-82 kev, the energy of the L photoelectron 
line will be 422-06 kev. The value given by Ellis was 421-98 kev but the last figures 
cannot of course be considered as significant. 
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The Elastic Behaviour of a Crystalline Aggregate * 


By R. HILU 
H. H. Wills Physical Laboratory, University of Bristol 


Communicated by N. F. Mott ; MS. received 17th December 1951 


ABSTRACT. "The connection between the elastic behaviour of an aggregate and a single 
crystal is considered, with special reference to the theories of Voigt, Reuss, and Huber and 
Schmid. The elastic limit under various stress systems is also considered ; in particular, 
it is shown that the tensile elastic limit of a face-centred aggregate cannot exceed two-thirds 
of the stress at which pronounced plastic distortion occurs. 


* This paper is a slightly extended version of a lecture delivered on 14th September 1951 at the 
Bristol Summer School on the Physics of Solids. 
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$1) INTRODUCTION 


s a method for calculating the elastic moduli of an aggregate of crystals, 
A Voigt (1928) proposed the averaging, over all possible lattice orientations, 
of the relations expressing the stress in a single crystal in terms of the 
given strain. Reuss (1929) proposed, on the other hand, the averaging of the 
relations expressing the strain in terms of the given stress. Voigt’s procedure 
contains the assumption that the strain is uniform throughout an aggregate, 
and Reuss’s that the stress is uniform. In the first model, the forces between 
grains could not usually be in equilibrium, while in the second the distorted 
grains could not fit together. In general, then, both derivations are approximate. 
It is found, however, that the measured moduli for the aggregate invariably lie 
between the Reuss and Voigt values, the latter being greater (Boas and Schmid 
1934). It is intended here to suggest an explanation, and to consider also other 
elastic properties of an aggregate. 


§2. RELATION BETWEEN CRYSTALLINE AND POLYCRYSTALLINE 
ELASTICITY 


Let the macroscopic stress and strain in the aggregate be denoted by S and E; 
for convenience of notation, these are regarded as vectors in a nine-dimensional 
space. Let o and € be the microscopic stress and strain at any point in any grain 
of the aggregate. Then, according to Bishop and Hill (1951a), if the aggregate 
contains a sufficiently large number of grains and is macroscopically homogeneous 


S.E={o.edV sete) 


integrated through a (large) unit of volume. 
In view of Maxwell’s reciprocal relation, the following identities are true at 
each point of the aggregate 


o.€+(o—o*) . (e—E)=o* .E+2(e—E) .a, 

o.€+(o—S). (e—e*)=S. e* +2(o—S).€. 
o* is the stress that would exist in a crystal with the local orientation and having 
the strain E; e* is the strain that would be produced in such a crystal by the 


stress S. Since the scalar product of any related stress and strain is essentially 
positive (being equal to twice the energy density) 


o.e€<o* .E+2(e—E).o 
and o.e€<S.e*+2(o—S)..«. 


It was shown by Bishop and Hill (1951 a, b) that the volume integral of the last 
term in each inequality vanishes. Hence, from (1) 


E. fo*dV, 
S.E< 
< eae Bi sta (2) 


The left-hand side of this equation is twice the actual energy density (as 
conventionally defined in terms of macroscopic quantities); the right-hand 
sides are twice the energy densities that would be calculated by the Voigt and 
Reuss theories respectively. 


The Elastic Behaviour of a Crystalline Aggregate 351 


For an aggregate which is macroscopically isotropic, (2) must reduce to 
the following symmetric (and automatically see coe expressions in the principal 
components of stress and strain 


a +E, + Es)? + aye —£,)< ae + Ey + Es)? +3 GyU(E; — E,)°, 
3 
K ag +S2+S3)*+ 56 aS, —S2)?< & (Sy + S2+S3)? + Emin bee 


K and G are the actual bulk and rigidity moduli, and suffixes V and R denote 
values calculated in the Voigt and Reuss theories from an evaluation of the 
integrals in (2). (Averaging the energy density, provided it is expressed in these 
respective forms, is equivalent to averaging the stress-strain relations for given 
.Stress or strain.) Since eqns. (3) are true for any stress and strain, it follows that 


Ka<K<Ky,; Gp<G<Gy. sieietanate (4) 
The corresponding estimates of Poisson’s ratio v and Young’s modulus E are. 
obtained by substituting the appropriate values of K and G in 


a pe Ge. ee iy Le 5 
ne 2} Ger Kee Go) ac SO ge OK ap dae e ©) 
From (4) it is evident that Ex, <E<Ey. 


It has thereby been proved that the Voigt moduli exceed the Reuss moduli, while 
the true values should lie between them. ‘The former conclusion is purely 
mathematical, but the latter is dependent on the truth of the statistical hypothesis 


embodied in (1). 
It was shown by Voigt that 
; 
OK y = (C41 + Coe + €33) + 2( C12 + Cog + C51), (6) 
15 Gy = (ey + Cog + C33) — (C12 + Cag + C31) + 3(Caa + C55 + G6) 5 i 
and by Reuss that 
1/Ky =(511 + S22 + 533) + 2(512 + 593 + $31); [ (7) 
15/Gp =4(551 + S22 + $33) — 4(S12 + S23 + 551) + 3(S44 + 555 + See). | 
It may be remarked that only nine of the 21 independent elastic constants for the 
single crystal occur in these formulae for the macroscopic moduli. 
In particular, for a body or face-centred cubic crystal 
Cy1 = C22 =C33, C2 =Cog= C31, 44 = "55 = Coe» 
with similar relations for the s coefficients, when the components of stress and 
strain are referred to the cubic axes. Also, 


Cy — yg =1/(S1- Siz), Ca F2Cy2=1/(Syy +2512), Cag=1/Sqy. +... (8) 
Equations (6) and (7) thus reduce to 
Ky =Ky=4(C1 + 2¢12), | (9) 


5 Gy =(€41 — C1g) + 3e4q, 5/Gp=4(S11 — S12) + 3544. i + Anat 


That K,=Ky is otherwise obvious since a hydrostatic stress produces only an 
isotropic volume change in face- or body-centred cubic crystals, and so uniform 
distributions of hydrostatic stress and strain in the aggregate are compatible.* 


* This would not be true of an aggregate of hexagonal crystals, for instance. 
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We may notice, in passing, that the Cauchy relation cy,=¢y, would make 
Gy =3K/5 and vy=}. Again, the condition for isotropy, ¢,1—¢)2:=2¢4,, makes 
Gp=Gy(=c¢y4); this is naturally to be expected since the aggregate is then 
microscopically homogeneous. Finally, for small amounts of anisotropy the 
difference between Gy and Gy is only a second-order quantity, since 
3 [2¢44— (C11 = C12)? 
Gy—Gp S (teu, wile eed ee (10) 


§3. NUMERICAL ILLUSTRATIONS 


Complete and reasonably accurate data are available for aluminium, copper, 
gold and «-iron. The c’s in table 1 are taken from Boas and Mackenzie (1950). 
Figures given by other authorities (e.g. Hearmon 1946) occasionally differ by 
about 1%. The tabulated values are for atmospheric pressure; according to 
Lazarus (1949) a superposed hydrostatic pressure of order 101° dyn/cm? increases 
the moduli by several per cent. 


Table 1 
Ci Cre Caa Gr Gy Krp=Ky &,R Ey VR Vy 
Al 1:08 0:622 0-284 0:26 0:26 0-78 0-71 0-71 0:349 0-348 
Cu 1-70 Lens 0-75 0-40 0:54 1-39 1:09 1-44 0:369 0-328 
Au 1:86 157 0:42 0-24 0°31 1-67 0-69 0:87 0-431 0-413 


a-Fe DEB ee aA Ll GO: 74 O39. 178 1-93 2:29 ¥ 0-313 0-280 
Unit=10" dyn/cm?. 


It will be observed that an aluminium crystal is much more nearly isotropic 
than the others. 

Table 2 shows the measured constants for the aggregates. Authorities differ 
appreciably, owing possibly to the presence of slight preferred orientations in the 
specimens, and to the inherent difficulty of determining v or K. The values 
adopted here have been chosen to conform to 1+v=£/2G. The agreement 
with (4) and (9) appears satisfactory, in view of the uncertainty in the data. The 
measured values of G fall almost midway between those of Reuss and Voigt. (We 
may contrast this with the conjecture by Bishop and Hill that the work done in 
plastic distortion may be calculated, to a good approximation, as if the microscopic 
strain were uniform.) Empirically, }(G_p+Gy) or 1/(GpGy) suggest themselves 
as good approximations. 


Table 2 
G K E v 
Al 0-265 0-74 0-71 0-34 
Cu 0-436 1°33 1-18 0-35 
Au 0-278 1-66 0-79 0-42 
a-Fe 0-808 1-59 2:07 0-285 


Unit= 10? dyn/cm?. 


Mention may be made of a theory proposed by Huber and Schmid (1934) 
and elaborated by Boas (1934, 1935). Their procedure for calculating 
Young’s modulus is tantamount to assuming that in any grain the stress is a 
pure tension parallel to some fixed direction, and of amount sufficient to produce 
a unit component extension in that direction. Effectively, a Young’s modulus 
is defined for each direction in the crystal and is averaged over all possible 
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‘directions. In the Reuss theory, on the other hand, it is the inverse of the modulus 
that is averaged. Evidently, therefore, the Huber-Schmid modulus must exceed 
the Reuss modulus. Calculations by Boas and Schmid (1934) for the four metals 
considered above show that the Huber-Schmid modulus exceeds Ex, by less 
than 0-3(Ey—E,). However, it has not yet been proved theoretically that the 
Huber—Schmid modulus is necessarily and invariably less than the true E or even 
less than Ey. 


ot. DHE PLASTIC wIM ra 


The elastic limit of an aggregate under a given kind of macroscopic stress 
(tension, shear, etc.) is taken here to refer to the stress at which local plastic 
deformation is initiated. The elastic limit is necessarily dependent on any 
locked-up stresses existing before the external load is applied. It seems probable 
that such initial stresses must tend to lower the elastic limit. Conversely, the 
reduction of such stresses by a mild annealing often serves to raise the elastic limit. 

In an aggregate of anisotropic crystals the elastic limit must also depend on 
the shapes of the grain boundaries, since these influence the local concentrations 
of stress. For this reason, we should expect a scatter in the measured elastic 
limits for macroscopically identical aggregates; other things being equal, the 
scatter should increase with the degree of elastic anisotropy in the constituent 
grains. 

It is possible to estimate an upper bound to the elastic limit, corresponding 
to what are presumably the most favourable conditions. Suppose that the 
Schmid law for plastic yielding applies; that is, yielding in a crystal begins when 
the component shear stress, acting over any slip plane and in any slip direction, 
reaches a critical value + (depending on the state of hardening). Suppose, 
further, that all grains are initially stress-free, equally hardened, elastically 
isotropic, and perfectly elastic until slip begins under the Schmid law. A 
uniform microscopic distribution of strain would then be produced by a uniformly 
applied external load not exceeding the elastic limit. Since all orientations are 
supposed to be present, local yielding begins in the grain(s) of ‘weakest’ 
orientation when the macroscopic maximum shear stress attains the value 7. 

Bishop and Hill (1951 b) have calculated that pronounced plastic distortion 
should first become possible in a face-centred aggregate under tension at the 
stress 3-067*; the writer has suggested the term ‘yield point’ for this stress. 
The tensile elastic limit, on the other hand, should be not greater than 27, that is, 
about two-thirds of the yield point. At loads between these two, the macroscopic 
strain (though partly plastic) is of elastic order (i.e. stress/F). Hence, until the 
yield point is reached, the load-extension curve deviates only slightly from 
the prolongation of the elastic line. ‘The corresponding ratio for pure shear 
is 3/5, very nearly. 

Extremely sensitive measurements are required for a reasonably close 
determination of the elastic limit, if the departure from the elastic line is due to 
slip in a few grains at most. Roberts, Carruthers and Averbach (1951) have been 
able to measure strains down to 2 x 10~®, and have reported elastic limits and 
yield points in several steels. For an annealed 0-17°%, carbon steel with only a 
slight drop in load following the upper yield point, they observed an elastic 
limit of 1-7 x 10° dyn/cm? and a lower yield point of 2-4 10° dyn/cm?. ‘The 

* Strictly, this value is for an aggregate whose elastic moduli are indefinitely large ; it should be a 
very good approximation for face-centred metals. 
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observed ratio was therefore a little over 2/3. This is in accord with the previous — 
remarks, since it appears justifiable to suppose that the ratio calculated for a 
face-centred aggregate is approximately true for mild steel. It was further © 
observed that the elastic limit had decreased to 1-5 x 10° dyn/cm? after 9°% strain, — 
while the yield point had increased to 3-5 x 10° dyn/cm?. ‘The ratio had thus 
fallen to 3/7, a circumstance which could be attributed to internal stresses 
introduced by cold-working. In the other steels, with various yield point drops, 
the ratio of elastic limit to lower yield point was between 2/3 and 4/5. 

Consider, now, the elastic limit under a pure hydrostatic stress. In an 
aggregate of face- or body-centred cubic crystals we have seen that the microscopic 
stress distribution is uniform in this case. Hence, if initially there are no locked-up 
stresses, the elastic limit under hydrostatic stress is indefinitely great, provided 
the Schmid law is valid (or, more generally, provided hydrostatic stress does 
not cause plastic distortion). On the other hand, in an aggregate of non-cubic 
crystals an external hydrostatic stress produces non-uniform microscopic 
stresses and strains, and so local yielding must eventually occur. Thus, there 
is an elastic limit in a non-cubic aggregate under hydrostatic stress, but, as Bishop 
and Hill have shown, no yield point within the range of validity of the Schmid 
law. Hence, both plastic and elastic components of the strain increase 
indefinitely as the hydrostatic stress is increased, while the rate of increase of 
strain remains of order & and decreases only gradually. 


§5. ELASTIC BEHAVIOUR DURING UNLOADING 


Finally, we may enquire for what range of stress an aggregate might be 
expected to unload elastically immediately after a permanent distortion. 
According to the theory proposed by Bishop and Hill (1951 b), the microscopic 
stresses in a face-centred aggregate during plastic distortion differ considerably 
from the macroscopic stress, so much so, that, on starting to unload, local slip 
on certain previously active slip directions would occur at once.* According to 
this model, then, no elastic range is to be expected in the first part of the 
hysteresis loop. In practice, an elastic range is usually observed, but the extent 
to which this is a reflection of the sensitivity of the apparatus seems not to have 
been investigated. 
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ABSTRACT. The eigenvalues and wave functions for the conduction electrons in metallic 
sodium are calculated for various points of high symmetry in the Brillouin zone. The 
method adopted is an extension of that used by von der Lage and Bethe in which boundary 
conditions are applied at a large number of points on the surface of the atomic polyhedron. 
Eigenvalues are obtained for states whose wave vectors lie at the ends of the 2-fold, 3-fold 
and 4-fold axes. The energy gap at the centre of a face of the zone is found to be 0°65 ev. 
A method of testing the results, of a more general nature than the empty lattice test of 
Shockley, is outlined and applied to certain cases. 


SiN TRODUCTION 

HE cellular method of determining wave functions and eigenvalues of 
electrons in crystals originated in the theory of Wigner and Seitz (1933) 
for the state of lowest energy of the conduction electrons in metallic 
sodium. Soon afterwards Slater (1934) developed the method to apply to states 
of higher energy, and he and his collaborators have made many calculations of 
the energy, as a function of the wave vector, for several metals. Unfortunately, 
the accuracy of these calculations is far less than that of the original Wigner—Seitz 
calculation. The reason is that the boundary conditions, which should be 
applied over the surface of the unit polyhedron, are much more complicated 
for states of higher energy than for the ground state, and consequently only 
approximate, and very inadequate, boundary conditions were, in fact, applied. 

A great improvement in the method was made by von der Lage and Bethe 
(1947, hereinafter referred to as LB), who showed that for certain special points 
in the Brillouin zone results of significant accuracy could be obtained. ‘This 
improvement was achieved by making use of the known symmetry properties 
of the wave functions at the particular points in the Brillouin zone considered. 
The theory of the determination of the symmetry properties of the wave functions 
was given in a fundamental paper by Bouchaert, Smoluchowski and Wigner 
(1936, hereinafter referred to as BSW). 

The notations used hitherto to specify an electronic state at a particular 
point in the Brillouin zone are not really satisfactory. LB deal only with the 
centre and six equivalent corner points of the body-centred cubic Brillouin zone 
and use Greek letters to denote states of specified symmetry. ‘This is. 
unsatisfactory because it would require a new set of letters for each point 
considered in the Brillouin zone: BSW use a capital letter to indicate the point 
in the zone and numerical subscripts to distinguish between the different states 
associated with that point. This is also unsatisfactory, because the symbol 
denoting a state gives no direct information about the symmetry of the wave 
function, and this is of interest in many physical applications, e.g. soft x-ray 
emission. We have therefore modified the BSW notation to remedy this defect ; 


the modification is described in §2. 
24-2 
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One purpose of the present paper is to extend the calculations of LB to — 
other points in the sodium Brillouin zone, in particular to obtain the energy 
gap at the centre of a zone face. We find this energy discontinuity to be 0-65 ev. — 
In §5 we obtain a correction to the eigenvalues, in the form of a surface integral ; 
apart from increasing the accuracy of the calculation of an eigenvalue, this 
correction can be used as a criterion of the success of any chosen set of 
approximate boundary conditions. Hitherto the only check on the accuracy 
of a calculation by the cellular method has been the empty lattice test introduced 
by Shockley (1937). This clearly is limited, in its application, to electronic 
states in real metals whose wave functions approximate, at least roughly, to free 
electron waves. The test described in §5, on the other hand, applies to any 
state however strongly the self-consistent field within the metal modulates the 
wave function. 


§2: THE SYMMETRY PROPERTIES OF THE WAVE LUNE TIONS 
The first step in the cellular method is to expand the wave function in a series 
of spherical harmonics, as follows: 


$= 2(A,, cos mé+ B,, sin mh)Pi” (cos O)Rife,7). ws eee (1) 
lm 
The radial wave functions R, are determined by the equation 
d?R, 2dR, (eee 
PE par aus + {e-V@)- 2 Ket pete cots (2) 


in which energies are measured in Rydberg units and lengths in units of the 
radius of the first Bohr hydrogen orbit. It is assumed in all this work that the 
potential energy V is spherically symmetrical within the unit polyhedron of 
the lattice. 

The wave functions associated with a particular point in the Brillouin zone 
fall into orthogonal sets, each set with a particular type of symmetry. If we 
select, therefore, a particular type, the number of coefficients in (1), up to some 
specified value of /, is in general reduced. The higher the symmetry of the type 
considered the smaller will be the number of coefficients remaining in (1) to be 
determined by the boundary conditions. 

The symmetry of a wave function, with specified wave vector k, must be 
stated with respect to the operations which transform the Brillouin zone into 
itself without changing k, or by transforming k into an equivalent vector. These 
operations form a group called by BSW the group of the wave vector k. Thus, 
to take an example, let k correspond to the point N in fig. (1), ie. a mid-point 
of a face of the Brillouin zone. The only equivalent k is that at the mid-point of 
the opposite face. ‘Che symmetry operations of this group of the wave vector k 
are therefore as follows: C,(z, Il) which denotes a rotation through 7 about the 
two-fold axis which is parallel to k; C,(7, L) which is a rotation through 7 about 
the two-fold axis at right angles to k (this transforms N into its equivalent point) ; 
C,(7), which is a rotation through 7 about the four-fold axis TH which is at 
right angles to k; J the inversion operator which reflects every point in the 
origin; and J times each of the above rotations. Including the identity 
operator £ there are, therefore, eight symmetry operations which will transform 
a Bloch wave function at the point N into another Bloch wave function belonging 
to the same k, 
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Since the Schrédinger equation is invariant with respect to these operations, 
it can be shown that the wave functions associated with the point N fall into 
eight types, each with a distinctive symmetry. In this paper we shall use a 
notation to indicate a particular symmetry type which is different from that 
used by BSW and which shows at once the number of true nodal planes through 
the centre of each atom by means of the well-known spectroscopic letters 
s, p, d, f, ete. We use the same capital letter as do BSW to indicate the point 
under consideration in the Brillouin zone (see fig. 1), but instead of a numerical 
subscript we use use the letters s, p, d, etc., as subscripts to denote the symmetry. 
Since in the crystal wave functions with the same number of nodal planes do not, 


Fig. 1. Brillouin zone for sodium; letters indicate points at which the wave functions 
have special symmetries. 


Table 1 
E (Oy: (zr) C, (7, II) Ce, (7, 1) J SC, (77) JIC, (77, 11) SC, (7, Ay) 
Ng 1 4 1 1 1 1 1 1 
N,? 1 == 1 —1 —4 1 —1 1 
N,? 1 1 —1 —1 —1 —{ 1 1 
N,? 1 4] —1 1 —1 1 1 —1 
Nat 1 =f 1 —1 1 —1 1 —1 
Na 1 = —1 1 1 —1 —1 1 
Na? 1 1 —1 —1 1 1 —1 =i] 
Ne 1 1 1 1 —1 —1 —1 a= 


in general, belong to the same type we have added a superscript (a number) to 
distinguish the wave functions with the same number of nodal planes. ‘Table 1, 
which is taken from the character table XII of BSW, shows the eight symmetry 
types belonging to the point N in the present notation. 

The number 1 in the table means that the wave function is symmetrical with 
respect to the operation shown, and --1 means that it is antisymmetrical. Thus 
the wave function :J(N,) has no nodal planes; the wave function /(N,!) a single 
plane node in the plane containing the axes C,(1) and C, (see fig. 2 (6)). These 
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are the functions with the two lowest energies at N for sodium; they correspond 
in the ‘extended k-space’ of elementary theory to the states on either side of 
the point N on the first Brillouin zone. It should be noted that the state N,?, 
which has a single nodal plane in the plane containing C,(1L) and C,(\j), 1s not 
degenerate with N,!. It is, in fact, a state of much higher energy which in the 
approximation of ‘nearly free electrons’ would have many additional nodal 
planes. However, in a real metal these additional planes are distorted, and only 
the node in the C,(1), C,(II) plane remains a true plane of antisymmetry, and 
this, therefore, defines the type. 

Table 2 gives the relation between the notations of BSW, LB and the one 
adopted in the present paper (HJ). This table, besides correlating the different 
notations, gives information which can only be obtained, after some calculation, 
from the character tables of BSW; for example, we see that the state denoted 
by BSW as I',, and by LB as y, is in the present notation [4!, i.e. a state whose 
wave function has two nodal planes through the atom. ‘This state, therefore, 
would behave rather like an atomic d state with regard, for instance, to the 
emission of X-rays. 


Table 2. Classification of Wave Functions in a Body-centred Cubic Crystal 


Equivalent points H, Mid-points of faces N Corner points P 
HJ BSW LB HJ BSW HJ BSW 
T's Ti o% Ng N; Pe Pi 
Tp Ths ) N,' NY’ PS PA 
Tq Tie ay Np? N,’ Ra Ps 
PGE IPs € N,° Ny’ Pe ie 
fs | BR B Na’ Ne Pj Py 

2 , 

y 25 Na’ N3 On trigonal axis 
g Ds’ : Na Ne A and equivalent 

Th Ti. Y Ng N2 type F 

; ype 

Tj ey B’ Ag Ay 
Ik Tr; oy Ap As 
Ag Ae 

Points on 4-fold axis A Points on 2-fold axis & Points on BZ face D, G 

HJ BSW HJ BSW HJ BSW 
Ag Ay Xs x, Ge G, 
Ay A; py 2 Gy! G3 
Aa A, x)” Da Gp? G,; 
Aa? Ay’ Xa mo Ga Gs 
As A 

1D% Dy 
Dy? Ds 
Da Ds 


In this paper we shall only be concerned with the states of relatively low 
energy. ‘These are : N,, N,'; H;, Hp, Ha’; P,, P,. Taking the known symmetry 
of the wave functions into account, as given by the character tables of BSW, 


the expansions for N, and N,! in spherical harmonics, up to 7<8, can be written 


#(N;) =Ry + ARP, (cos 0) + BR,P, (cos 0)+ CRgP, (cos 6), ...... (3) 
¥(N,') = RP, (cos 0) + AR;P3 (cos 0) + BRP; (cos 0) + CR,P, (cos 6), 
AON fe (4) 


where A, B and C are, in each case, unspecified constants. 
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For any particular type of wave function we are at liberty to choose the polar 
axis, with respect to which @ is defined, in the most convenient way in order to 
make the expansion as simple as possible. For the two states N, and N,' we 
take as the polar axis C,(II) (cf. fig. 2 (b)), and the dihedral angle ¢ is measured with 
respect to the plane C,C,(II). 

The expansion (3) was obtained as follows: From table 1 we deduce that 
the wave function #(N,) is symmetrical with respect to reflection in each of the 
coordinate planes C,(L) C;(II), CxC2(II) and C,C,(1). Thus, with the polar axis 
as chosen, we see that there can be no term in the expansion in sin md. Similarly 


Fig. 2. The atomic polyhedron for sodium (see text). 


there can be no odd spherical harmonic, i.e. P," (cos 0) cos mf in which either 
m or lis odd. However, besides the terms shown in (3), others, of the above 
form with even m, could occur, e.g. P,? (cos 9) cos 24. These terms are not 
excluded by the known symmetry properties but have been omitted, in our 
approximate expansion up to 1<8, because we use the empty lattice test to 
estimate the accuracy of our energies, and the terms in question do not occur in 
the expansion of the corresponding free electron wave function. This 
emphasizes the weakness of the empty lattice test and shows the need of a 
criterion such as that given in §5. 

For the wave functions associated with the point H (fig. 1) we take C, as the 
polar axis and measure ¢ relative to one of the planes containing two fourfold 
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axes (see fig. 2(c)). From the character table I of BSW, or from the data in LB, 
we deduce the following forms for the functions belonging to the states 
eh He; Hy. 


~(H,) = Ry + AR,{P, (cos 9) + 745P,' (cos 9) cos 44} 


+ BR,{P, (cos 9) — 335P 64 (cos 8) cos 4h}, we eee (5) 
¥(H,) =R,P, (cos 6) + AR3P3 (cos 0) + BR;P; (cos @) + CR,P, (cos @), ....+. (6) 
ala a = R, P,* (cos @) cos 26+ AR,P,? (cos 6) cos 2¢ 

+ BR,{P,? (cos 9) cos 26 + x45 P66 hes 0) c0s:6}.9 1". ae (7) 


In the expansion for the wave function 4(H,) we have omitted the harmonics 
P;4 (cos 0) cos 44 and P,4 (cos @) cos 4¢ for the same reason that similar terms 
were omitted from %(N,) and %(N,1). The state H, is triply degenerate; the 
particular wave function given by (6) is the one in which the nodal planes are 
parallel to the plane containing C,(L) and C,(I}). 

The group of operators belonging to the point P differs from those already 
considered (the groups of N and H) in that it does not contain the inversion 
operator J. Hence the wave functions associated with P do not fall into sets of 
odd or even parity. We shall write 

oP) Sue, and==d(P )=u oy | eee (8) 
where in both cases u and v are complex functions which have even and odd 
parity respectively. 

The group of P contains 24 operators, and if we add the products of the 
inversion J and the operators of P we obtain the full cubic group of 48 operators. 
u, is invariant with respect to each of these 48 operators and has therefore the © 
symmetry of I... From table XI of BSW we can read off the characters of P, 
(i.e. BSW, P,), and if J is represented by —1 for the term of odd parity v, the 
characters of the remaining 24 operations will be just those of this table (line 1) 
multiplied by —1. By comparing this set of characters for 48 operators with 
table I of BSW we see that v, has the symmetry of type I" (LB, 8). Ina similar 
way we can show that u, has the symmetry of I, and v, the symmetry of I). 
From this knowledge of the symmetries of the wave functions we can determine 
the expansions in spherical harmonics, and we find 


u,=R,y + AR,{P, (cos 8) + 4gP,4 (cos 4) cos 44} 
+ BR,{P¢ (cos 9) — 335P,* (cos 8) cos 64}, Ss... ss es (9) 
vs = CRP? (cos @) sin 26 + DR,{ P;? (cos @) sin 24 + 74,P,8 (cos 0) sin 6¢}, 
esa (10) 
and for the P,, type 
u, = CR, P,? (cos @) sin 24 + DR, P,? (cos 8) sin 2d 
+ FR,P,? (cos 6) sin 24+ GP,§ (cos @) sin 6d, ...... (11) 
vp = RP, (cos 0) + AR3Pz (cos 0) + BR;P; (cos 0) + K P;4 (cos 0) cos 44. 
sagan (12) 
In the expansions (9), (10), (11) and (12) the angles 6 and ¢ are measured with 
respect to the same axes as were used for the H types. 
From considerations of symmetry alone no further information about the 
wave functions can be obtained than is given in the expansions of this section. 


The unknown constants along with the eigenvalues must now be determined by 
the boundary conditions. 
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§3. THE EXACT BOUNDARY CONDITIONS 
Let n denote one of the fourteen vectors which may be written {1, 1, 1} and 
{2, 0,0} with respect to the fourfold axes C,. A face of the atomic polyhedron 
may then be denoted by the vector }an which is the displacement from the 
origin to the centre of the face. If R, denotes the displacement along the 
normal from any point on the face —}an to the face fan, then clearly R, = 4an. 
In accordance with Bloch’s theorem we have, therefore, 


o(r+R,)=exp(zk.R,)b(r), = 2...0. (13) 
n. Vi(r+R,,) = —exp (#k.R,)n.Ve(r). °  ...... (14) 


Let X, denote a rotation through 7 about the normal to the face n through the 
centre, then 


Ae ta IC a) ee ee (15) 


where J is the inversion operator. For the points N and H in the Brillouin zone 
the operator J is contained in the group of the wave vector and, therefore, every 
wave function associated with these points has either odd or even parity, 
i.e. Ji=+y. For such points in the Brillouin zone the boundary conditions 
may, therefore, be written 


Me Bar exp (kraRi juin) aV Linen’ ais oe: (16) 
n. VX ob(r)= 4 exp (kk Rijn Viir) eS (17) 


where the upper sign is taken in both equations when ¢ is of even parity and the 
lower sign when ¢ is of odd parity. 

The boundary conditions (16) and (17) are easier, both to visualize and to 
apply, than (13) and (14) because they now relate to points on the same face 
of the atomic polyhedron. We calculate, as an example, the boundary conditions 
-over the faces shown in fig. 2(d) for the function %(N,). The coordinates of the 
point N in the wave vector space, with respect to the axes C, shown in fig. 1, 
are (z/a, w/a, 0) and, for the hexagonal face with centre A (fig. 2(d)), 
R,,=(4a, 4a, 4a) referred to the axes C,. Thus k.R,=7 and, since #(N,) 
has even parity, we find X,4=—%. ‘Thus, referring to fig. 2(b), we have 
ux = —ig, where the subscript letters denote the points at which the function 
is evaluated, and since the function is symmetrical with respect to a reflection in 
the line CZ we have 4,=-—y,. Hence the wave function must be anti- 
symmetrical with respect to AQ and, therefore, vanishes along this line. 
Evidently the normal derivative is symmetrical with respect to AQ. 

The hexagonal face with centre D corresponds to n=(1, —1, 1), with respect 
to the C, axes, therefore k.R,=0. Hence the wave function is symmetrical 
with respect to DQ and the normal derivative antisymmetrical. ‘The face MQL 
is given by n=(2,0,0), therefore k.R,=7 and thus X,4=— 7 so that the 
function is antisymmetrical with respect to BQ, which is a nodal line. ‘The 
normal derivative is symmetrical with respect to BQ. 

For the face FNC we have R,=(0,0,a) and, therefore, k.R,=0. ‘Thus 
X,w=% and X,’=—', where the dash denotes differentiation with respect 
to the normal to the plane. In this case X,, the rotation through 7 about C,, 
is a member of the group of the point N, and we know that for the N, type the 
wave function is symmetrical with respect to this rotation. ‘Thus the first 
condition is automatically satisfied and the second necessitates ;’=0 over the 
whole face. 
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In this way the exact boundary conditions, for the types considered at 
N and H, have been obtained and are shown in table 3. #’ denotes the normal 
derivative of % to the face which contains the line specified. 


‘Table 3 
Fig. 2 (6) Line AQ Line DQ Line BQ Face FNC 
op p op ip ip p op ip 
N, a s s a a s 0 
Ny} s a s a s a 0 
Fig. 2 (c) Line AQ Face QBL Face NCX 
p p p p ib p 
Hy s a 0 0 
Hg? s a 0 Q.* 
Fig. 2 (d) Line AN | Face NCX 
p yp op yr 
1S a s 0 


s and a denote symmetry and antisymmetry with respect to the line indicated; 
0 implies the vanishing of ¢ or %’ over the whole face. 

The boundary conditions for the types P, and P,, require more detailed 
consideration. Consider P, first. For the hexagonal face n=(1, 1,1) we find 
k.R,=37/2 and thus JX,%= —2. Now let u,=a,+78, and v,=A,+7z,, where 
as, Bs A, and pu, are all veal functions. Thus, since Ju,=u, and Jv,= —v,, we 


find easily X wep, and KS E L eeee (18) 


Now «, stands for the expansion (9) in which the unknown constants are the 
real parts of A and B and uw, for the expansion (10), in which the constants are 
the imaginary parts of C, D. It will now be seen that the two equations (18) 
give exactly the same equation to determine the eigenvalue, by elimination of 
the unknown constants, and we need, therefore, only consider one of them. Thus, 
without loss of generality, we put 8,=A,=0 and write J(P,)=u,+7v, where 
u, and v, are now real functions and may be taken to be the expansions (9) and (10), 
in which all the constants are real numbers. 

Since u, has the symmetry of [,, for which the wave function transforms 
into itself under each of the 48 operations of the full group, it is only necessary 
to consider the region OARNCX of fig. 2(d), which is one forty-eighth of the 
whole polyhedron. The boundary condition (18) now implies that u,(T) =2,(T’) 
where, as shown in fig. 2(d), [and ‘I’ are diametrically opposite points on che 
hexagonal face. Since the plane OVA is a symmetry plane for v, (which is the 
LB f type) we can write u,(T) =v,(T”), which we may express = saying that 
u.—v, 18 antisymmetrical with respect to the line AN. The full boundary 
conditions are summarized as follows: 


‘Table + 
Fig. 2 (d) Line AN Face NXC 


(us—Vs) (ts — 5.) Us Ve! 


os 


a (P.) a Ss 6) 9) 
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The symmetry of the P,, function may be defined with respect to fig. 2(c). 
For u, the plane NQB is a nodal plane, i.e. a plane of reflection antisymmetry, 
whilst the planes XCAZ and OBLZ are planes of reflection symmetry. For 
v, the planes XCAZ and XNQB are planes of reflection symmetry whilst the 
plane OBLZ is a nodal plane. The boundary conditions determined in the 
same way as for P, are given in table 5. 


Table 5 
Fig. 2 (c) Line AQ Face NXC Face QBL 
(up—Up) (Up’—vy’) Uy Up. Up Up’ 
us (Py) a s 0 0 0) 


As in table 3, primes indicate differentiation along the normals to the planes 
indicated. 


§4. APPROXIMATE BOUNDARY CONDITIONS AND 
NUMERICAL RESULTS 

From the exact boundary conditions, given in the previous section, we can 
obtain a set of particular relations which may be used to determine the constants, 
in the expansions of the wave functions, and the eigenvalues. Since we take 
only a finite number of terms in each expansion, there is inevitably an 
arbitrariness in the choice of the particular boundary conditions to be applied. 
For each type we have tested many sets of difterent boundary conditions. Those 
which give the best results in the empty lattice test are considered to be most 
suitable for the determination of the energies in the real metal, at least in the 
case of metallic sodium. 

For the wave function %(N,) table 3 shows that the following equations must 
hold (subscript letters refer to fig. 2 (b)): 


(1) fot%z=9, (it) oo =z, (itt) Hg = 9, (iv) py +ym=9, 
(v) by’ =m (vi) py =9, (vil) x'=0, (vill) py =0, 
(ix) pq =0, (x) $a’ =9, (xi) bp =0. 
Where the letter refers to a corner point of the polyhedron the particular plane 
to which the normal derivative applies can be seen at once by reference to table 3 
and fig. 2(b); for example the derivative in y’ =0 must be with respect to the 
normal to the plane NXC, whilst (x) must refer to the plane YDFQ. 

The procedure followed in applying these boundary conditions was as 
follows: A value of the energy was chosen and the radial wave functions 
determined. In the empty lattice test the radial wave functions are Bessel 
functions and their values can be taken from existing tables. Any one of the 
above eleven conditions leads to a linear equation in the three unknown 
constants 4, B and C of expansion (3). If we select three of the conditions, 
e.g. (ix), (i), (v), we can determine A, B and C for a given value of the 
eigenvalue «. Hence we can calculate, for example, the numerical value of yp’ 
of condition (xi). ‘This calculation is repeated for a few values of ¢ from which, 
either numerically or graphically, the value oi e« can be found which makes 
%p vanish. We have theretore obtained the wave function, i.e. the constants 
of (3), and its eigenvalue by satisfying four of the boundary conditions. 
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Table 6 gives the eigenvalue of the N, state as determined by different sets of © 
boundary conditions, and also the percentage error in the empty lattice test. 
In the ‘empty lattice’ the eigenvalue is exactly 27?/a?. A percentage error of © 
(0-036, for example, means that the eigenvalue, calculated by the expansion (3) — 
with approximate boundary conditions, is 19-732/a” in place of 19-739/a’. 


Table 6 
Conditions used Error in E.L.T. (%) er oe Me 
(ix), (x), (xi), 0-036 —0-319 
(ix), (x), (xi), (ii1) 0-036 —0-318 
(ix), (x), (xi), (iv) 0-046 —0-318 
(ix), (x), (xi), (v1) 0-046 —0-315 
(ix), (i), (viii), (xi) 0-030 2595315 
(xi), (ait), (vi), (1) 0-091 —0-+299 
(xi), (iii), (vi), (1x) 0-036 —0-310 
(xi), (iii), (vi), (x) 0-032 —0-299 


Adopted mean value for the eigenvalue of Ns=—0-315 Ryd. 


In a similar manner a set of conditions, satisfied by the wave function %(N,1) 
at various points on the octant of the polyhedron shown in fig. 2 (), can be deduced 
from table 3. The conditions we have selected as most suitable are given below: 


(i) Yo=Yin (ii) oo +x! =0, (iil) ’ =0, (iv) Yo= hn 
(v) py +¢m' =9, (vi) op =9, (vii) (O/Or)g=0, (viii) (0/28)g =9, 
(1x) py’ =0. 


Choosing sets of four of these conditions, the results given in table 7 were obtained 
for the eigenvalues of the N,]! state. 


Table 7 
Conditions used Error in E.L.T. (%) penn a a 
(ii), (vii), (viii), (vi) 0-286 —0:327 
Gi), Cae G@id. ds) 0-0456 —0-270 
(vii), (x), Gv), Gi) 0-0456 — 0-269 
(vii), (ix), (iv), (v) 0-0912 —0-267 
(ii); Gx), (v); (vii) 0-0760 —0-268 
(ii), (ix), (v), (@) 0:1267 —0-267 
(ai), (ix), (v), (iv) 0-0862 —0-267 


Adopted mean value for the eigenvalue of Np'= —0-268 Ryd. 


The first line of table 7 shows that a poor empty lattice test leads to an eigenvalue 
appreciably different from the rest. It has been included in order to serve as an 
example for the correction procedure which is developed in the next section. 
The eigenvalues of five other states: H,, FH, Hat; Ppt and) Pichavey been 
calculated in the same way. In each case several sets of boundary conditions 
were applied, those leading to small errors in the empty lattice test being 
considered most reliable. It is unnecessary to repeat the details of these 
calculations, and we give only our results in table 8. 
The eigenvalue for T, has been taken from LB and was included in table 8 to 
complete the list of accurately known eigenvalues for sodium. 
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For any particular state the differences of the eigenvalues in the empty lattice 
test, resulting from the application of different combinations of boundary 
conditions, are very rarely above 1%, and in the majority of cases less than 0-1%, 
as exemplified by the two tables 6 and 7. Plots of the approximate wave 


‘Table 8 
Biss Eigenvalue for Na Free electron value relative 
(Ryd.) to [T's (Ryd.) 

ie —0-608 
Ns —0-315 —0-303 
Np' —0-268 

p —0-0135 —0-:012 
Ha? +0-0172 
He +0:0997 
ie —0-110 —0-152 
Pa +0-020 


Energy (electron volts) 


H 
i (Q) Mang oS (b) (c) 


Fig. 3. Energies as functions of the wave vector : (a) [110] direction, (0) [111] direction and 
(c) [002] direction. Detailed figures show the deviations from the free electron energies 
(broken lines). (c) also shows a comparison with the energies at H as determined by LB. 
Energies are measured in electron volts relative to the energy of the I's state. 


functions in the empty lattice test along various selected directions showed 
very close agreement with the known true standing waves of a free electron. 
It may be noticed, for comparison, that in the original work of Slater and his 
collaborators errors in the empty lattice test of up to 20°, occurred frequently, 
as was shown by Shockley. 
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It is clear that the accuracy of the results obtained for any real metal must 
depend upon the accuracy with which the radial wave functions R,(e, 7) can be 
obtained by numerical integration. Since our result for H,' does not agree 
with that obtained by LB, who give — 0-036 Ryd., a check was made in our radial 
wave functions by using them together with the boundary conditions employed 
by LB. Exact agreement was obtained with LB’s numerical result. The 
difference, therefore, arises from the different boundary conditions which were 
applied. A plot of the function, along the line NL (fig. 2(c)), obtained by the 
LB boundary conditions, showed much less symmetry about the line AQ than 
the function obtained by our point boundary conditions. Strong evidence will 
be given in the next section for believing that our value is the correct one. 

The quantities which are of most interest physically are usually differences 
between the eigenvalues, e.g. the energy gap across a face of the Brillouin zone. 
It is for this reason that high accuracy in the eigenvalues is essential if the results 
are to have a practical significance. Consider, for example, the two states 
N,and N,!. The adopted means for the eigenvalues are — 0-315 and — 0-268 Ryd. 
respectively, with the greatest error in either case about + 0-002, i.e. less than 1%. 
The energy gap, however, is only 0-047 Ryd. and, in the least favourable case, the 
error might be +0-004, which is 8%. Expressed in electron volts the energy 
gap at the centre of a face of the Brillouin zone is 0-65ev. It is, perhaps, 
noteworthy that Butcher (1951), by comparing the absorption power of sodium 
with the calculated inner photoelectric absorption, obtained a value of 
0-323 ev for the Fourier coefficient of the potential. In the approximation of 

‘nearly free electrons’ it is well known that this Fourier coefficient is equal to 
one half of the energy gap. 

Figure 3 shows graphically the results of the arekoae calculations, and also 
the way the various states link together within the Brillouin zone. In the present 
example, sodium, the way the states join together is fairly obvious; in general 
it must be deduced from the compatibility relations given by BSW. 


§5. CORRECTION TO THE ENERGY BY MEANS On 2Ai 
SURFACE INTEGRAL 

In the present application of the cellular method to the case of metallic 
sodium we have used the criterion of the empty lattice test to estimate the 
accuracy of the true eigenvalues in the metal. In this section we show how a 
correction to the energy may be obtained, and how this may be used to provide 
a test of the accuracy of any energy determined by approximate boundary 
conditions quite independently of the empty lattice test. 

Let ¥%(r) be the approximate wave function which satisfies the equation 


Vebotlegoe Vito =O 9 DSL aewES (19) 
exactly, but fits the boundary conditions only approximately; and let (x) be 
the function which satisfies the boundary conditions exactly and the equation 


Ve a (e227) pa eal oe id aaah (20) 


We now calculate the difference (e—«)). Multiply (19) by % and (20) by %&; 
subtract and integrate over the cell. Then 


(€ — €o) orb dr = [LV hq — oV hb) dr = SUV itbo —YoVrib) dS... (21) 
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by Green’s theorem, where V,, denotes differentiation along the normal to dS. 
The surface integral of (21) can be reduced further when we know the exact 
boundary conditions which % satisfies. We consider the case of H,! in detail; 
others can be treated similarly. In this case the surface integral of (21) is 16 times 
the value taken over the three planes shown in fig. 2(c). Consider the surface 
NCZLQ. From table 3 we see that %(H,') is symmetrical and V,:/(H,') 
antisymmetrical about the line AQ. Let any point K on the upper half NCAQ 
be denoted by r and let r’ be the reflection of K in AQ, i.e. K’.. We now define 
two functions of position on NCAQ as follows: 


Play rn ee aes (22) 

ets Vaud tat VUE ea ames. (23) 
Thus the boundary conditions can now be expressed as f(r)=g(r)=0 at every 
point of the surface NCAQ. Let f((r) and g,(r) be defined in the same way 
as f and g except that the approximate wave function ys, is used in place of x. 


If, now, S, denotes the area NCAQ and S, the area AZLQ, then, by using (22) 
and (23), we may write 


I gs PVntho— YoY) 2S = | (PV nitho — PoV nh) dS 
+ | , {leo —Vaslo) + (ho—Fo)Vail} 4S... (24) 
= | .Ws0—foVub) aS, 


where we have made use of the fact that # is exactly symmetrical and V,% exactly 
antisymmetrical with respect to AQ. 
Over the surface QBL, which we denote by S3, we find V,% =0, from table 3. 


Hence 
[PV nbo—YoV) dS = | bVitodS, sae (25) 


and similarly the contribution from the surface NCX, which we denote by Sy, 
may be seen to be Jfg,bVy~)dS. Adding the contributions from the three 
surfaces together we obtain from (21) 


Heo) [Yo dr= | Wgo—foVat) dS+ | bVitodS+ | YVapodS. ....(26) 


This equation is exact ; but we now approximate by replacing & by % throughout 
(26). This amounts to neglecting products such as (%—y%)g) over S,, and 
(us —)Vyx% over Sz and S,. These are all quantities of the second order, since 
£, on S, and V,y% on S, and S, tend to zero as ys tends to %. Thus, if zs is 
normalized to unity, we have as a first order approximation 


B(e—€o) = |, Yoko—foVaho) AS + |, PoVaodS, sees (27) 


where the integrands are now known functions. 

Equation (27) may be used as a test of the accuracy of any eigenvalue <p, 
or wave function y, obtained by approximate boundary conditions, and is quite 
independent of the use of the empty lattice test. We give below two numerical 
examples of the use of this equation. 

The first entry in table 7 gives —0-327 Ryd. for the eigenvalue of N,,’, which 
is quite exceptionally far from the mean value. It can be shown that the 
correction for the state N,' is given by’ 


He-eo)= | ,(foso—FoVato) 25+ | toVatodS, evens. (28) 
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where o, denotes the surfaces NCAQ+FNQD+QBL (fig. 2(6)) and o2 the 
surface FXCN. 

We have, therefore, calculated the integrals of (28) using the corresponding 
wave function and found that «—e,=0-043, giving a corrected eigenvalue of 
— (284, which is much nearer to the mean value. The difference between this 
corrected value and the mean is probably due to computational errors in the 
evaluations. of the integrals; no attempt was made to obtain these with great 
precision. 

As a second example we consider the difference between the eigenvalue of Hy" 
as calculated by LB and by ourselves. Using the wave function obtained by 
applying the ‘average’ boundary conditions of LB which, it was verified, 
correspond exactly to the eigenvalue given (—0-036), the integrals of (27) were 
evaluated. It was found that « —«,) =0-060, from which we deduce « = + 0-024 Ryd. 
The value obtained by our boundary conditions was 0-017 (cf. table 8). The 
difference between the corrected LB value and ours, viz. 0-007 Ryd., is again 
probably due, in the main, to lack of precision in the estimation of the surface 
integrals. 
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Nuclear Spins of Samarium 147 and 149 


Recently Murakawa and Ross (1951) reported the nuclear spins of samarium 147 and 149 
as being 5/2 in each case. ‘They based this on their observation of a maximum of six 
hyperfine structure components in the optical spectra of the separated isotopes. By means 
of paramagnetic resonance measurements on the ethyl sulphate of naturally occurring 
samarium we find instead that the spins are 7/2. The experimental evidence is presented 
below. 

Bleaney, Elliott and Scovil (1951) showed that the spin-lattice relaxation time of 
samarium ethyl sulphate was too short at 14° kK for the resonance spectrum to be observed. 
We have extended these resonance experiments to 4° K, where one electronic transition is 
found. 

The spectrum is completely resolved in the concentrated ethyl sulphate and is essentially 
the same in a crystal diluted with lanthanum ethyl sulphate in the ratio Sm:La=1:100. It 
consists of a single strong line flanked by sixteen satellites. The strong line is due to the 
even isotopes 144, 148, 150, 152 and 154 which have spin zero. The sixteen satellites fall 
into two groups of eight as shown in the figure. The two groups are identified as belonging 
to the isotopes 147 and 149, whose relative abundances are 15-07+ 0-15°% and 13-84+ 0°14% 
respectively (Inghram, Hayden and Hess 1948). Since the hyperfine structure should 
consist of (2/-+1) lines from each isotope, the spin is shown to be 7/2 in each case. 'The 
intensity of the 147 lines relative to the strong line is 2:-69+ 0-1%%, which agrees favourably 
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with a relative intensity of 2:65°%, predicted by the known abundance and a spin of 7/2. The 
relative intensity of the hyperfine structure of isotopes 147 to 149 is 1-:10+0-02, which 
also agrees with the abundance ratio. Thus both spins are undoubtedly 7/2. Since this is 
in disagreement with the conclusions of Murakawa and Ross, we suggest that they did not 
observe transitions between levels with sufficiently high J to show the full multiplicity of the 
hyperfine structure. 

Our measurements show that the ground electronic state of the samarium trivalent ion 
in this salt is a Kramers’ doublet whose spectroscopic splitting factor g (in the dilute salt) is 
0-596 + 0-002 parallel to the hexagonal axis and 0-604+ 0-002 perpendicular to the axis. 


149 i l ! | | | | | 
147 


a 


‘The hyperfine structure of Sm 147 and 149 arising from the M= +4~>M= —} electronic transition 
in samarium ethyl sulphate. In this oscillogram, which was taken at an angle of 12-6° to the 
axis, the marked asymmetry of the spectrum results mainly from the large second-order effect 
near the axis and partly from the non-linearity of the sweep. 


The hyperfine structure is described by the usual spin Hamiltonian : 
Bye Sa SP gg 11) AS BO el, yl tee ac s1(I+ 1)} . 
‘The values of the constants are (dilute salt) 


Sm 147 A=0-0060+ 0-0001 cm-1 B=0-0251+0-0001 cm7} 
Sm 149 A=0-:0049+ 0:0001 cm- B=0-0205 + 0-0001 cm7 


An upper limit of 4x 10~* cm can be assigned to the P of either isotope. The precision of 
the measurements is shown by the following comparison of the observed and calculated 
positions (in gauss) of the eight hyperfine lines of isotope 147 in the perpendicular direction 
at a wavelength of 3-279 cm. 


Observed 13833 12819 11856 10930 10054 9228 8368 7598 
Calculated 13835 12830 11865 10933 10043 9190 8379 7607 


This agreement could only be obtained with J=7/2. The same holds true in the parallel 
direction and for the other isotope. 

The ratio of the nuclear magnetic moments, given by the ratio of the widths of the hyper- 
fine groups, is 4147/u149=1-:222+ 0-008, in fair agreement with the value of 1:177+0-015 
given by Murakawa and Ross. If this discrepancy is genuine it may be due to the fact that 
we are dealing exclusively with a 4f electron configuration, whereas they were observing 
transitions between levels involving s electrons. A nuclear size effect may be expected in 
the latter but not in the former case. It is interesting that the addition of two neutrons 
(N=85 and 87) produces a considerable decrease in the magnetic moments of the nuclei, 
as in the case of Nd 143 and 145 (VN=83 and 85) (Bleaney and Scovil 1950), where it was even 
more marked. 

The actual values of the nuclear magnetic moments and upper limits for the nuclear 
electric quadrupole moments have been calculated by Elliott and Stevens (1952), and are 
given by them in an accompanying letter. 


The Clarendon Laboratory, G.S 
Oxford. H. E. D. Scovit. 
15th February 1952. 
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Exiott, R. J., and Stevens, K. W. H., 1952, Proc. Phys. Soc. A, 65, 370 

IncurRaM, M. G., Hayden, R. J., and Hess, D. C., Jr., 1948, Phys. Rev., 73, 180. 
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The Paramagnetism of Samarium Ethyl Sulphate: Theory 


In an earlier letter (Elliott and Stevens 1951) we reviewed the paramagnetic resonance 
observations made on the rare earth ethyl sulphates. A resonance had not then been 


a ee 


found in the samarium salt because of its low g-values and short relaxation time. In ~ 


the accompanying letter Bogle and Scovil (1952) report new observations on this salt and it 


appears that the lowest doublet is, as predicted, approximately |*H;;.; J,=+1/2). The - 


next multiplet is near 1000 cm-! and small admixtures from this considerably affect the 
g-values. In fields of Dz, symmetry the lowest Kramers’ doublet will be approximately 


cos 6 |J=5/2, J,=+1/2)+'sin @ |J=7/2, J,= + 1/2) 


where 0 is related to the crystal fields. Then 


a= 5 {ooe 6+ 2 sin* 6+-34/10 sin 2 ak 


3 5 
f= ={3 cos? 6 — sin? 6 — xv 10 sin 2 af 
and the ratio of the hyperfine structure constants is 


3 

{6 cos? 6 +- _ sin? @—34/10 sin 2 at —s 2 cos? 6 -+- sin’ 6— 10 sin 2 at 

148 34/10 76 3 Pe, : 
{18 cos? §— yr sin? 6 ae sin 2 a} 7 ‘6 cos? 9+ TE sin? 9-+ 77105" 2 a} 
where 5=13/45. The hyperfine structure arises from the interaction of the nuclear 
magnetic moment with the orbital and spin magnetic fields of the electrons. A characteristic 
feature of the rare earths is that the spin contributions (here given by the terms in 8) are 
small. With 0=4° the comparative values of the experimental and theoretical values of 
yp g, and A/B are 


byl As 


Si 81 A/B 
Theory 0-66 0-66 0-24 
Experiment 0-595 0-605 0-24 


The agreement is fairly good and can be improved by including the smaller admixtures 
from otherJ levels. Even when this is done there remains a discrepancy for what seem the 
most reasonable crystal fields. The details of this theory will be reported in forthcoming 
papers on the whole group of salts. 

The theory gives p (1/r?)=-+1-52 x 10°B, where 1/r® is averaged over 4f electron wave 
functions and expressed in A~*. (1/r*) is difficult to estimate, but from the spin—orbit 
coupling constant and the value for Eu given by Sternheimer (1950), a reasonable value 
would be 40 to 50. With the values of B given by Bogle and Scovil, we find for 14’7Sm, 
|| =0-68 + 0-1, and for Sm, |u| =0°55+0-1 nuclear magnetons. The nuclear quadrupole 
moment Q as defined by Mack (1950) is given by Q(1/7?)=+7:2 x 104 P, with P in cm—. 
The upper limits of |P| =4 x 10~4 set by Bogle and Scovil give upper values for |Q| of 0-72 
in 147Sm and 14°Sm. 

The values of the magnetic moments are about 50% larger than those obtained by 
Murakawa and Ross (1951), even when 7/2 is substituted for the value of 5/2 which they give 
for the nuclear spin. 

Spins of 7/2 are more consistent than 5/2 with the nuclear shell model (Mayer 1950). 
“Nd has 83 neutrons—one more than a magic number—and since it has spin 7/2 (Bleaney 
and Scovil 1950) it is apparently the nuclear 2f,;, state, not the 1hg;2, which is preferred at 
the beginning of this shell. ‘The corresponding proton case, ?°°Bi, has spin9/2. So far other 
isotopes (°Nd, *“’Sm, “°Sm) with three and five neutrons outside a closed shell have spin 
7/2. Moreover *’Er which has the 7/2 and 9/2 shells full except for one neutron, also has 
spin 7/2 (Bleaney and Scovil 1951). The addition of more neutrons fills these shells and 
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1Yb and '8Yb have spins of 1/2 and 5/2 (Mack 1950). Thus one might tentatively 
predict that all isotopes with an odd number of neutrons between 83 and 99 will have spins 
of 7/2, or just possibly 9/2. 


The Clarendon Laboratory, Ree heer. 
Oxford. K. W. H. STEvENs. 


15th February 1952. 
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Proton Magnetic Resonance in Solid Cyclohexane 


The nuclear magnetic resonance absorption spectrum and the spin-lattice relaxation 
time J, have been measured for the protons in polycrystalline cyclohexane between 100° k 
and its freezing point (279-6° Kk). The second moment (mean square line width) of the 
measured spectrum is shown as a function of temperature in fig. 1. Measured values of T, 
are shown in fig. 2. 

Below 150° k the spectrum does not change, indicating that the molecular lattice is 
effectively rigid. 'The mean experimental value of the second moment between 100° kK and 
150° K is 26-0+0°5 gauss”. A theoretical value of 27:3 + 1-5 gauss” has been obtained using 
the formula of Van Vleck (1948). The agreement between measured and calculated values 


Spin-Lattice Relaxation Time 7; (sec) 


Second Moment (gauss?) 
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Fig. i. Fig. 2. 


provides general confirmation of the molecular configuration assumed in the calculation, 
namely D3q4 symmetry with tetrahedral bond angles of 109° 28’, C-C bond lengths of 1:54 & 
and C-H bond lengths of 1:10 A. The main uncertainty in the calculation arises from the 
lack of exact knowledge of the crystal structure below 186° k, at which temperature there 
is a first order phase transition. 

Between 155° k and 180° k the absorption line narrows, the second moment falling to 
6-4 gauss’. ‘The molecular motion most likely to cause this narrowing is the re-orientation of 
each molecule about its triad axis. Calculation of the effect of such motion on both the 
intramolecular and the intermolecular contributions to the second moment lead to a 
theoretical value of 6:1+1°0 gauss? in satisfactory agreement with the experimental value. 


25-2 
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If the re-orientation process can be described by a single correlation time 7, the relation 
between 7, and 7, should be of the form derived by Bloembergen, Purcell and Pound 


fee 1/T,=Cy[rei(1 +0" 70°) + 27,/(15-40°r.) 1, ee (1) 
where w/2z is the radio frequency in use (23:4 Mc/s). Below 186° k w?7,?>1 since Ty, 
increases rapidly with decrease of temperature; relation (1) thus approximates to 
T,=2w"7,/3C,. Taking te=7T» exp (E/RT), a plot of log T, against 1/T yields for E, the 
height of the barrier restricting re-orientation about the triad axis, the value 11 + 1 kcal/mole. 

The line width-temperature curve has been reconstructed in the manner described by 
Gutowsky and Pake (1950), and can be fitted using the same value of E. This enables the 
value of C, to be found, namely 3 x 10° sec~®. An approximate value of C, may be calcu- 
lated by assuming that the relaxation of each proton is predominantly caused by re-orienta- 
tion relative to the other proton in the same methylene group. The relation C,=3yh?/10b® 
due to Bloembergen et al. (1948), which strictly applies to isotropic re-orientation as in 
liquids, may then be used. (y is the nuclear gyromagnetic ratio of the proton and 3 is the 
interproton distance.) A value of 5 x10°% sec~? is obtained in satisfactory agreement with 
the experimental value. 

Below 160° k JT, is shorter than is to be expected from this re-orientation mechanism 
alone. ‘Thus, below 160° k, this relaxation process ceases to be dominant, other processes 
having comparable or greater efficacy. 

The solid—solid phase transition at 186° K is accompanied by discontinuous changes in 
both second moment and relaxation time. The reduction of second moment to 1-4 gauss? 
implies that the molecules are now able to re-orient about other axes besides the triad, and 
since JT, is increasing with temperature the frequency of re-orientation 1/277, must now 
exceed the radio frequency 23:4 Mc/s. Thermal data indicate that there must be a relatively 
large increase in molar volume at the transition ; this in turn must reduce the barriers hinder- 
ing molecular re-orientation, and thus accounts for the greater molecular freedom observed. 

The decrease of T, above 240° k is presumably associated with the further increase in 
molecular freedom which causes the final reduction of second moment between 220° k 
and 240° k to less than 0-01 gauss?. ‘The absence of a measurable discontinuity in either 
the second moment or in JT, at the melting point emphasizes the considerable molecular 
freedom of the molecules just below the melting point. 


Department of Natural Philosophy, E. R. ANDREW. 
The University, ; R. G. EabEs. 
St. Andrews, Scotland. 
18th February 1952. 
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Wavelength Resolution of X-Ray Proportional Counters 


In a recent letter on a “ Comparison of Geiger and Proportional Counters for Intensity 
Measurements in X-ray Diffraction”, Arndt and Riley (1952) state that the proportional 
counter’s “‘ response to any particular characteristic wavelength employed can be made as 
sharp as required”’. It should be pointed out, however, that the wavelength resolution of 
proportional counters is fundamentally limited by statistical fluctuations in the number of 
primary ion pairs formed and in the amount of gas amplification. On the other hand, 
Arndt and Riley were unable to discriminate completely against harmonic radiation according 
to the account given of their experiments, whereas the writer has previously shown (Lang 
~ 1951) that separation of CuKa radiation from its first harmonic can be complete. 

The purpose of this note is to deal quantitatively with the case, important in practice, of 
detecting CuK« radiation in the presence of slightly softer fluorescent radiations which can 
be excited by the CuKa rays. 


‘Gade 
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‘The variation of pulse height from CuKa radiation was found to follow closely a gamma 
distribution function. The experimental points and the fitted curve are shown in fig. 1. 
Distribution curves for the other radiations can be computed from the CuKa distribution 
since the variance of the distribution is proportional to the energy of the radiation. The 
ratio of variance to the square of the mean pulse height obeys the relation o?/p?=C/N 
where N is the mean number of ion pairs formed before gas multiplication and C is a constant. 
For CuKa radiation N= 298, taking the value of the energy required per ion pair formation 
in argon to be 27 ev (Valentine 1952). With this value of N the constant C is found to 
equal 1-1. The distribution curves shown for CrKa, MnKa, FeKa and CoKa radiations 
are all normalized. The FeKf curve has an area 0-13 of the Kacurves. The experimentally 
observed subsidiary peaks centred 3 kev below the main peaks, which are due to loss of 


FRACTION OF PULSES COUNTED. 


PULSE HEIGHT ENERGY EQUIVALENT IN keV. DISCRIMINATOR SETTING ENERGY EQUIVALENT IN KeV. 
1. Pulse height distribution curves for argon Fig. 2. Computed bias curves for argon proportional 
proportional counter. counter. 


argon K fluorescent radiation in the absorption process, are not shown here since they 
contain only about 4°% of the number of pulses in the main peak (Lang 1951). This fraction 
will perforce be lost when discrimination against slightly softer radiation is attempted. 

From measurements on these curves, or on bias curves computed from them (fig. 2), the 
degree of discrimination obtainable is readily calculated. It will be seen that CrKa, 
MnkKaand FeKa can be almost completely separated from CuKa radiation. Discrimination 
at 7-5 kev will pass 85°%, of the CuKa radiation main peak and only 10% of the CoKa 
radiation. FeKf radiation can be largely excluded but not so CoKf radiation (only the 
FeKf pulse distribution is shown in the figure). 

The help of Mr. D. R. Cox and Mr. D. A. East of the Cambridge University Statistical 
Laboratory in computing the distributions is gratefully acknowledged. 


Crystallographic Laboratory, A. R. Lane. 
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20th February 1952. 
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The Ultra-Violet Spectrum of a Diatomic Hydride 


In the course of a search for the emission spectrum of KF, which should consist of 
fluctuation bands in the ultra-violet region, we have recently examined the emission from 
a Schiiler discharge-tube, running in hydrogen, in which the heated hollow cathode contained 
KF. No KF bands were observed, but a many-line spectrum in the region 2150-2650 A was 
sometimes developed with some intensity. [he spectrum appeared to be the same as that 


374 Letters to the Editor 


obtained under similar conditions by Woods (1943), who suggested that it might be the 
25 —2]] system of HF+. No analysis of this system seems yet to have been published. On 
the assumption that the emitter was HF+, we photographed the spectrum when KHF, was 
used in place of KF in the cathode. The same many-line spectrum appears strongly in this 
source, but it is somewhat capricious of excitation, and a good deal of difficulty was 


experienced in getting plates adequate for a preliminary rotational analysis. However, ° 


satisfactory spectrograms on a Hilger large quartz instrument have now been taken. The 
results of the analysis are unexpected. 

Most of the strong lines, about 230 in all, have been assigned as P or R branch lines to one 
of five bands constituting a v’=constant progression. ‘The bands are degraded strongly to 
longer wavelengths, and no heads are apparent. The values of the constants derived in the 
preliminary analysis are given below. 


Ue Vo Bie Dis 

0, 0 46380 11-920 1-47 x 10-3 
Ors 44345°5 11-261 159 

Oe 424112 10:524 iePKy, 

0, 3 40779°5 9-665 to 

OVA: 392913 8-718 Dae 


F’ (J)=4:016, J(J +1)—1:89, x 10-4 J2(J +1)?+-4-.x10-° J3J(+1)?. 
o =484/D)i7=-1170;em* 

Be. =12:215 056,03) 0-047 474-3) - 

Gy" =2174-(0" + 4) —64-0,(0" + 9)?—3-69(0" +. 
In the lower rotational states the terms in HJ*(J-+1)* become of increasing importance as 
v” increases : because of the magnitude of these terms the values of D” given above are 
only approximate. [he vibrational quantum numbers are assigned on the assumption that 
the band at shortest wavelength is the 0, 0 band: it may prove later that these values of 
v’, v” should be increased. 

Further work designed to identify the emitting molecule is in progress. The constants 
do not fit with those of any known state of a diatomic hydride, but the number of possibilities 
is not large. The value of B,” is such that the emitter must almost certainly be the hydride 
of a first-row element (see, for example, Douglas and Herzberg 1942). Although other 
possibilities cannot yet be eliminated, it is probable that the transition is 1X—1X, for no 
Q branches have been found, which indicates that AA =0, and no spin or A-doubling effects 
have been observed up to J=30, although lines about 1 cm™ apart are clearly resolved. 
If the transition is +U—14, HF* is ruled out, and of the possible emitters remaining, CHt+ 
or OH? are the most likely. 


Physical Chemical Laboratory, A. D. CauntT. 
University of Oxford. R. F. Barrow. 
28th February 1952. 
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Derivation of the Velocity of Second Sound from 
Maxwell’s Equation of Transfer 


Starting from the Boltzmann equation for the distribution function, Ward and Wilks 
(1952) have recently derived the relation cy=c//3 for the velocity of second sound in a gas 
of phonons moving with velocity c. In this note we use an argument of a similar type to 
calculate c, for a general energy spectrum E=E (p) (where p is the momentum of an excita- 
tion), and compare the resultant theory with that of Landau (1941). 

The most convenient starting point is Maxwell’s general equation for the transfer of any 
quantity Q per unit volume. ‘This reads (cf. Jeans 1940, § 191) 


a— i a ee 
ve + div (VQ) + iV. grad: O= Aoistong O55 | eistsseie (1) 


Re 
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where V is the velocity of a particle or excitation carrying QO, V its acceleration, and the bars 
denote values integrated over the distribution function f. For instance, 


O=55|{[ Orp* sin # dp aba = jp | Ofer. bcd seal (2) 


We restrict our attention to the case of plane waves propagated along the x direction. 
Assuming that the only forces acting on the excitations are those exerted during 
collisions, V=0. If Q refers to either energy or momentum, Agojisions Q=0, so that 


We cannot assume that the distribution function f is equal to its undisturbed symmetrical 
value fy>=1/{€ exp (E/RT) $1}, because the right-hand sides of (3) and (4) would then 
vanish while the left-hand sides would not. Instead, we assume that there is a small 
distortion of the distribution function, and therefore put f=f)+/,, where f, is the unsym- 
metrical part. Since v and p are parallel vectors, averaging over angles gives p,,v,, f>4pufo, 
while v,,f+v,fi>vf’ and p,,f>p,fi—pf’, f’ being a measure of the distortion of the distribution 
function. Thus (3) and (4) become 


i) fi] 1a () 
—_ 2 = Se if’ p2 pet 3 eee 7 13 
= | Efe dp = | Bef'p GIDE G, duct 0 (5) nee | veo dp =| rp GprMnte (6) 


The elimination of the unknown quantity f’ is particularly simple for a phonon gas, for 
then E=pc and v=c, and the integrals on the left-hand sides of (5) and (6) are identical, 
while the integral on the right-hand side of (5) is just c? times that on the right-hand side of 
(6). Thus @/@t2=4c? 0?/dx?, leading to a velocity cp=c/1/3. 

For a general energy spectrum, however, it is both simpler and more instructive 
to compare the equations with those of Landau’s theory, rather than attempt to eliminate f’ 
directly. The integrals on the left-hand sides of (5) and (6) are ordinary thermodynamic 
functions, for 


4 
7 | Bho dp=U, B=) nike Weve hs (7) 


the internal energy per unit volume, while putting v= 0E/ ép (cf. Dingle 1951, § 4, 1952, § 5), 
it is easily shown that, provided E (p)>oo as p> , 


4 
= a | ofa” Be Pah Ae Vi I (8) 


the free energy per unit volume.* For small temperature fluctuations, and neglecting 
thermal expansion, 
dU/at=TaS/et; oF /ex=—S dT] ox, 


so that (5) and (6) become 


os 0 4Ar te eee 
Ee ee ’p2 —S—=——|fpdp. ....-. 10 
ag ade as | Bef’ FM antes 3 (9) ae amie p a) 
These equations may now be compared with the following arising in Landau’s theory : 
as C) OT 0 f pnpvn 
eet Seem yo ew oe Te —S eaees tae sommes. ee pre tele) 0 Ye 12 
at aa On) oe ax = Pa (12) 
In (11) and (12) wv, is the ‘drift velocity’ of the excitations in the x direction. Since by (8) 
oF 4m fo 5 Aare {ous 
S=-— = —— J v= =———— | hvu pdp aca a aee 13 
anand sale ae? a cart one oe, 


the last step following from the fact that fo=fy (E/T), agreement between our equation (9) 
and Landau’s equation of entropy conservation (11) is obtained only if 


Pa=ahaPeEs 8 eee ease (14) 


* If the particles or excitations are conserved, the left-hand side of (8) is equal to F'— NE), where 
N is their number and E, their degeneracy energy. The argument is unchanged. 
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Equation (14) shows that the literal interpretation of v, as ‘drift velocity’ is not really 
necessary. Instead, it may just be regarded as a parameter representing the degree of 
distortion of the distribution function. 

Introducing into the right-hand side of (10) Landau’s effective mass-density in the form 
(cf. Dingle 1951, eqn. (4.7), 1952, (79)) 

4x ( , So 4a _ 

Pu=— 355 p oR {re Gp) | Ree (15) 
it is seen that complete agreement between (10) and (12) is achieved only in the limit 
Ps->P, i-€. Pp=p—ps—>0, which requires T--0. Put another way, our equations (9) and (10) 
would lead to a propagation velocity given by c.?= T:S?/Cp,, where C is the specific heat per 
unit volume, as compared with Landau’s relation c= TS?p,/Cppy obtained from (11) and 
(12). The explanation of this discrepancy is probably that any description in terms of the 
excitations alone, such as that presented here, is adequate only if their effective mass is 
negligible compared with the inertial mass of the total assembly. Apart from this defect, 
the present treatment based on Maxwell’s equation of transfer is, so far as it goes, essentially 
equivalent to Landau’s more complete theory. 


Royal Society Mond Laboratory, R. B. DINGLE. 
Cambridge. 
14th February 1952. 
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Diffusion in Gaseous Helium at Low Temperatures 


Progress in recent years in the separation of "He and measurement of it sproperties have 
increased interest in theoretical estimates of these properties. Much has already been done 
in this direction in the calculation of critical constants and vapour pressure (de Boer and 
Lunbeck 1948, Lunbeck 1950), the equation of state (de Boer, van Kranendonk and Compaan 
1950), and viscosity (Buckingham and Temperley 1950, de Boer and Cohen 1951) at low 
temperatures, where quantal effects are most important. This note adds information about 
the diffusion coefficient and thermal diffusion ratio in ?>He—*He mixtures at 5°K and below. 

The calculations use mainly an interatomic potential of the form introduced by 
Buckingham and Corner (1947) in dealing with the equation of state at higher temperatures. 
This is V(r)=— ef(o), o=r/ro where €, 79 are respectively the depth and position of the 
minimum of l’(r), and the shape function is given by 


fo)=fio (1+Bo) —fye—“e—0), 


where a, B are parameters, and f,=«a/{«(14-8)—6—8B}, fo=—1+(1+)f,. Strictly, the 
above form is used only for o greater than 1; when o<1 an exponential factor is introduced 
which has the effect of eliminating gradually the terms in o~®, a8. Although the effect of 
varying « and 8 would be interesting, the values «=13-5, B=0-2, known to be reasonably 
good, have been adhered to in the present work. It has also been assumed that 
€=14-04 x10 erg, r9=2:943 A, these values being based partly on a comparison 
with observed second virial coefficients of “He at low temperatures (Buckingham, Hamilton 
and Massey 1941). De Boer and his collaborators have used a potential of the usual 
Lennard-Jones type in their work. 


The table contains the following results : 

(1) The coefficient of diffusion, D,., for a mixture of He and ‘Heat a pressure of 1 mm Hg. 
The values are small, but for temperatures above 3°K they are about twice those predicted 
by classical theory (Hirschfelder, Bird and Spotz 1948). (2) The thermal diffusion ratio 
kp for a mixture containing a molar fraction y, of "He. For potential (1) and y,=0°3, 0-5, 
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0-7 results are given in columns (a), and it will be seen that the variation with y, is asym- 
metric, being more heavily weighted on the side y, >0-5. For comparison, some results 
obtained using quantal theory for rigid spheres of diameter 2-1A are given in column (6); 
also, the classical value for rigid spheres, which is independent of temperature. (3) The 
theoretical values of the viscosities, .(*He) and .(4He), which were used in the calculations of 
kp. ‘The values for *He were derived from potential (1), and supersede the preliminary 
results of Buckingham and 'Temperley (1950); those for “He are based on a slightly different 
potential used in earlier work, though the difference is unlikely to be important. 

All these represent first order approximations only; second order corrections, though 


possibly increasing as T becomes very small, and largest in kp, should not exceed a few per 
cent. 


—kp x 10? 
° D T (He) (4He) 

T (°x) ae (a) (a) (a) (6) F oetee 

(cm?/sec) 03 0-5 0-7 0:5 (10- c.g.s.) 
0-25 Ono) —9-8 —10-0 2235 0-52 
0-5 0-028 —0:52 —0-68 —0-64 2:86 3-62 1E9O: 
O75 0-059 +1:25 +1-66 +1-56 4-90 3225 
1-0 0-098 2°16 2°85 2:68 PNG 6:25 3-48 
1-25 0-144 229 3-03 2°87 foo3 3-61 
1:5 0-198 199 2-65 2-52 1-98 8-69 3-84 
2-0 0-325 1:88 2:49 2-36 1-83 10:3 4-82 
2 0-478 1:78 2°34 2:20 17 ite 6°05 
30 0-655 Oo) 2:20 2-03 1-65 11:8 7°47 
3-5 0-855 1-63 2:09 1-90 12-1 9-08 
4-0 1-07 1-60 2-03 1-83 12-5 10-8 
4°5 1232 1-64 2:06 1-83 13-0 22 
5-0 1-58 1-80 2:22 1-93 soy 13-0 


For rigid spheres (classical theory), —kp X 10?=3-09, y,=0°5. 


It is perhaps worth noting the sudden change in sign of kp near 0:5°K. The general 
behaviour of kp, as might be expected, departs considerably from classical theory at these 
low temperatures. Thus the calculations of Hirschfelder, Bird and Spotz, using the Lennard- 
Jones potential, indicate that kp is positive over a temperature range from about 4°K to 
9-5° (assuming e/k ~10°k), and attains a positive maximum of 0:0020 around 6°xK. 

Full details of this work will be published later. The calculation of transport properties 
at higher temperatures is also well advanced. 


University College, London. R. A. BucKINGHAM. 
10th March 1952. R. A. SCRIVEN. 
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Intensity Distribution of Bremsstrahlung from Beta-Rays 


The jemission of beta-rays is accompanied by internal bremsstrahlung, external 
bremsstrahlung and characteristic x-rays. This note is concerned only with the external 
bremsstrahlung. The intensity distribution, as a function of quantum energy, has been 
calculated approximately from the theory given by Heitler (1944), and the results have been 
checked by means of absorption measurements on the bremsstrahlung from **P. 

The intensity distribution of external bremsstrahlung for an electron which has a par- 
ticular kinetic energy E has been given by Heitler (1944, fig. 14). For moderate energies 
the shape of the intensity distribution is roughly independent both of the atomic number of 
the material through which the electron is passing and of the kinetic energy of the electron, 
and can be represented by the equation 


I (hv)=4—3hv/E, (0<hv<E, Ebeing constant), ...... (1) 
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The intensity of bremsstrahlung at a particular quantum energy hy, radiated from the whole 
track of the electron, may be obtained by integrating eqn. (1) with respect to E from the 
initial kinetic energy E, down to hv; this gives 


ih) =By44 (\-#) 3 Pho et Jie (2) 
Eo v 


The distribution of eqn. (2) is plotted as curve A in fig. 1, where the function I(A), (A=Av/Eo), 
has been normalized so that the area under the graph, which is the total intensity of brems- 
strahlung radiated from the whole track of the electron, is unity. Since the shape of the 
distribution given by eqn. (2) is independent of the value of Ey, curve A applies to electrons 
of any initial kinetic energy up to a few Mev. ‘To find the intensity distribution of brems- 
strahlung from a radioactive isotope, eqn. (2) must be integrated over the beta-ray spectrum. 
This has been done numerically for *?P, using the shape of the beta-ray spectrum given by 
Siegbahn (1946). The resulting intensity is plotted as curve B of fig. 1, where in this case 
A\=hv/Eg and Ep is the maximum beta-ray energy. Curve B has also been plotted so that 
the area under-the graph, which is the total intensity of bremsstrahlung radiated by the 
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Fig. 1. Intensity distribution of brems- Fig. 2. Absorption in lead of bremsstrahlung from °2P. 
strahlung from a single electron (curve A: The curve is calculated from theory, the dots are 
A=hyv/E,) ; and from a beta-ray spectrum measured values. : 

(curve B: A=hy/Es). 


isotope, is unity. Since the shape of curve B is independent of the value of Eg, up to a few 
Mev, the curve gives the intensity distribution of bremsstrahlung from any isotope whose 
beta-ray spectrum has the same shape as ®2P. ' 

In order to check these calculations, measurements were made of the absorption in lead 
of the bremsstrahlung from **P. The geometry of the experiment was good. The source of 
radiation consisted of 10 g of red phosphorus (B.D.H. ground), irradiated in the pile at 
Harwell, and contained in a thin-walled brass box. The total radiation from such a source 
consists of external bremsstrahlung, internal bremsstrahlung and characteristic X-rays, plus 
any gamma-radiation due to contamination. The total intensity and the intensity distribution 
of the internal bremsstrahlung has been calculated by Knipp and Uhlenbeck (1936), and 
measured by Wu (1941) and by Mandanski and Rasetti (1951). For the source used an this 
experiment the total intensity of the internal bremsstrahlung is about 20° of the intensity 
of the external bremsstrahlung, and the intensity distribution is not very different from that 
of curve B. Hence the internal bremsstrahlung will have a negligible effect on the shape 
of the absorption curve. 'The characteristic x-rays are so soft that they will be almost 
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completely absorbed in the wall of the counter. The intensity of gamma-radiation due to 
contamination was calculated from a chemical analysis of the phosphorus and should have 
been negligible. However, the first absorption measurements revealed the presence of a 
strong gamma-ray of about 1:75 Mev quantum energy. Fortunately the contaminating 
gamma-radiation decayed with a half-life of about 15 hours (it was probably due to 24Na), so 
that a week later its effect was negligible. 'The absorption curve in lead was now measured 
again, with the results shown by the dots in fig. 2. The probable error was 2°, in most cases. 

The curve of fig. 2 was calculated from the intensity distribution of fig. 1, curve B, 
taking into account the variation of the counter sensitivity with quantum energy (Renard 
1948, down to 0-1 Mev, Saurer 1950, below 0-1 Mev) and the initial filtration of the brass 
box. The agreement is good beyond 0-5 mm of lead, but the measured points indicate a 
rather greater proportion of very soft radiation than was calculated. 


Radiotherapy Department, S. J. Wyarp. 
University College Hospital, 
London W.C.1. 
18th February 1952. 


HEITLER, W., 1944, The Quantum Theory of Radiation, 2nd edn. (Oxford: University Press). 
Knipp, J. K., and UHLENBECK, G. E., 1936, Physica, 3, 425. 

MAnpanskI, L., and Rasetti, F., 1951, Phys. Rev., 83, 187. 

RENARD, G. A., 1948, 7. Phys. Radium, 9, 212. 

SaAuRER, H., 1950, Helv. Phys. Acta, 23, 381. 

SIEGBAHN, K., 1946, Phys. Rev., 70, 127. 

Wu, CHIEN-SHIUNG, 1941, Phys. Rev., 59, 481. 


REVIEWS OF BOOKS 


An Introduction to Experimental Physics, by H. G. JERRaRD and D. B. McNEILL. 
Pp.576. (London: Harrap, 1951.) 25s. 


Anyone who has had the responsibility of organizing an elementary physical laboratory 
is aware of two important problems: (i) Are students to be provided with manuscript 
directions appropriate to each experiment and the apparatus which is available for its 
performance ? Or are they to be expected to purchase a textbook of practical physics, 
the exercises in which approximate fairly closely to what can be done with the available 
apparatus ? (ii) Bearing in mind the limited number of sets of apparatus for a particular 
exercise, how are the experiments performed by the students to be linked with the theoretical 
treatment of the subject in the class-room or lecture theatre ? 

The book under review undoubtedly provides an exhaustive collection of experimental 
exercises suitable for an intermediate course in physics; but for this very reason, in size and 
price it compares with several well-known textbooks of theoretical intermediate physics. Ata 
price of twenty-five shillings it is certainly value for money by modern standards; yet it 
is doubtful whether all the students using a particular laboratory could be expected to 
possess a copy, in addition to equally expensive and more necessary books which cover the 
theoretical aspect of the subject. It does not seem, therefore, to be a substitute for “ manu- 
scripts’. Moreover, few students would cover adequately all of its 180 experiments in 
the time usually available for practical work. Inevitably, therefore, selection is necessary ; 
and the book appears to be more suitable as a volume of reference in the laboratory than as 
a manual to be possessed by every student. 

The problem of linking a student’s experimental and theoretical work is met by giving 
a fairly complete preliminary account of the principles underlying the group of related 
experiments which form the contents of each particular chapter. ‘Thus the volume is 
self-contained and on the whole orderly. It seems strange, however, that a fairly 
full treatment of pendulums, including the compound bar pendulum, in Chapter VIII, 
should precede the verification of the elementary laws of dynamics by means of Fletcher’s 
trolley in Chapter XIII. This again suggests that the book is most suited to be a work of 
‘reference. 
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Criticism of the book must be made on a number of points of detail. Though the 
diagrams are clearly drawn, a few are inaccurate or misleading. For instance, fig. 30 
suggests an extremely orderly packing of the molecules of a liquid. Surely, unequal, 
randomly oriented forces could be shown in equilibrium. Fig. 31 shows two distinct 
normals at a single point on a curved liquid surface! The treatment of the telescope 
(pp. 331-5) suffers in that the diagrams show no rays, other than paraxial rays, passing 
right through the system. Fig,300, illustrating the forces on the needle of a dip circle, in 
two planes at right angles, is obviously in error. The points O, A, C, B should appear in 
perspective as the corners of a parallelogram. There is also a certain amount of unnecessary 
duplication of diagrams. For example, 267 repeats 246; 268 repeats 247; 302 repeats 281. 

In the section on magnetism and electricity the reviewer feels that a more critical attitude 
might well have been adopted towards a number of well-known but not very valwable 
exercises. Hibbert’s magnetic balance appears to be somewhat out of place in a modern 
textbook, particularly without any discussion of the mutual demagnetization of the magnets ; 
and, in any case, the apparatus would be unworkable if the suspended magnet were ‘ supported 
at its centre on a knife-edge’, as stated (p. 421), or worse still as it is shown in fig. 284. 
The location of neutral points is discussed without reference to Owen’s method except 
in a short footnote (pp. 410-4). "The moments of magnets are compared by the ‘sum and 
difference ’ method, again, without reference to the precaution which is essential in order 
to avoid demagnetization (p. 428). Ohm ’s law is ‘ verified’ by voltmeter and ammeter, 
and only a footnote suggests that the experiment is worthless except as a very elementary 
illustration of the law (p. 460). The ‘sum and difference’ method of comparing the 
e.m.f.’s of cells is given without any reference to the well-known objection to this exercise 
(p. 481). The Wheatstone net is given in every case without a reversing key (pp. 493-507). 
Presumably it is only by an oversight that we read on p. 467: ‘‘ The absolute value of the 
volt is 108 ergs”’; while the section in which this statement occurs gives, to say the least, a 
misleading impression of the present status of the ‘ international’ volt, ampere and ohm. 

However, the book is clearly written, well-printed, and stoutly bound. Full use is 
made of graphical methods of representation, including the plotting of logarithms. The 
inclusion of ‘ worked examples’ sets a high standard for the student’s notebook, which the 
reviewer hopes will be copied. The defects that are noted above can easily be remedied 
in a future edition, and the book should find a place in many physical laboratories in 


schools and colleges, where it will serve a useful purpose as a singularly full and complete 
textbook of experimental physics. GSD: Xe 


Differential Equations, by H. B. Puituirs. Pp. viiit+ 149. 3rd Edition revised. 
(New York: Wiley and Sons; London: Chapman and Hall, 1951.) 24s. 


This is a most attractive introduction to the subject, intended mainly for physicists and 
engineers. While providing some formal instruction in the solution of simple ordinary 
differential equations, most emphasis is placed on the setting-up and solution of differential 
equations arising in Applied Mathematics, Physics, Chemistry, etc., such physical principles 
as are required being briefly described in the text. 

Chapter I deals with those first-order equations which can be solved by separation of 
the variables, with a wide variety of applications. 

Chapter II completes the work on first-order equations, exact, linear and homogeneous 
equations being discussed. ‘The integrating factor for the general linear equation is merely 
stated and verified. A direct derivation would take very little extra space. The alternative 
method of solving the same type of equation given on p. 44 so obviously comes to the 
same thing that it could well be omitted. The chapter ends with a brief description of 
simultaneous equations and solution by change of variable. 

The simpler types of second-order equation, i.e. those in which x or y is explicitly 
absent and homogeneous equations; are dealt with in Chapter III. Theoretical discussion 
is confined to eight pages only, but there are many worked examples, including a section 
on the motion of dynamical systems of varying mass. 

The final chapter treats linear equations of the nth order with constant coefficients. 
in an elementary way, operational methods receiving only a brief mention. 

In addition to about seventy illustrative examples worked out in full in the text, there: 
are nearly five hundred exercises, for which answers are provided at the end of the book. 


A. N. GORDON... 
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The Interpretation of X-ray Diffraction Photographs, by N. F. M. Henry, 
H. Lipson and W. A. Wooster. Pp. ix+258. (London: Macmillan, 
2951 ) ees. 


The appearance of this book marks a stage in the development of x-ray crystallography. 
A considerable number of books have been written on the subject, but this is the first text- 
book written specifically for the student pursuing a systematic course of study. X-ray 
crystallography has become an essential tool in many types of research, both fundamental and 
applied, and is increasingly used in industry. The need for the training of x-ray crystallo- 
graphers has been recognized in recent years and steps taken to provide it. This book will be 
of very great assistance in this task. 

The book covers the interpretation of x-ray diffraction effects from crystalline material, 
but does not deal with space group theory or the derivation of atomic positions from a list of 
corrected intensities. It has an excellent system of cross reference and a large number of 
original papers are referred to. The photographs and charts are produced to scale so that 
the student can carry out each step of the interpretation himself. A number of very useful 
tables are given in appendices. 

it is sufficient recommendation to say that this book is an essential textbook for any 
serious student of the subject, and, since it has no competitor, the remainder of this review 
can safely be devoted to suggestions for improvements in the next edition, which is certain 
to be required very soon. 

In the first place there is a certain lack of balance which is possibly due to the book’s 
triple authorship. The interpretation of Laue photographs is dealt with in great detail 
and occupies 16 pages, together with four more in the chapter on crystal setting. On the 
other hand, the interpretation of Weissenberg photographs is covered in 14 pages, and the 
measurement and correction of intensities (including a discussion of structure factor theory) 
also receives only 14 pages. In the latter case, the geometrical factor for equi-inclination 
Weissenberg photographs is not considered at all, and, in general, practical details of the 
methods of correction are dealt with only by references to original papers. In other places 
the text appears to be condensed to a degree which makes comprehension unnecessarily 
difficult. In view of the fact that the indexing of Laue photographs is only very rarely 
required, it is difficult to justify the relative space given to this subject. 

The second main criticism concerns the interpretation of oscillation photographs. In 
the process of setting a crystal, the recognition of salient features of a tilted reciprocal lattice 
may be a major step. ‘This recognition is facilitated if one is in the habit of viewing the 
reciprocal lattice and sphere of reflection with the x-ray beam coming towards the observer. 
In the book under review, however, these are presented with the x-ray beam going away from 
the observer. 

The problems of determining the azimuthal orientation of a crystal and the B-angle of a 
monoclinic lattice are insufficiently developed. A section on the design and use of camera 
arcs could well have been included in Chapter 4. 

There are two points of presentation which the reviewer would like to see incorporated. 
‘The effect of the Laue conditions for diffraction can be most easily visualized by the use of 
three intersecting cones whose axes coincide with the crystal lattice axes. It seems a pity that 
this should have been left out. 

The second point concerns the use of the reciprocal lattice in interpreting effects due to 
the ‘ white’ radiation. ‘The method adopted is to keep the reciprocal lattice fixed and vary 
the size of the sphere of reflection. However, for some purposes it is desirable to keep the 
sphere of reflection fixed in size and associate a weighted radial line with each characteristic 
point of the reciprocal lattice. This line is limited towards the origin by the short wave- 
length cut-off, and fades out on the other side of the characteristic point. ‘The production of 
Laue streaks on constant p curves, and the measurement of the end positions of a streak in an 
oscillation photograph on (90°—¢) curves, can be easily visualized in this way. 

In the chapter on fibre photographs the problem of interpreting a complicated fibre 
diagram, such as is given by some minerals, is not really tackled at all. The starting point 
for interpretation is given as the ‘ indexing of the power lines’. In the case of mineral fibres, 
the occurrence of a fibre photograph is usually the only thing that makes indexing possible at 
all, and the problem is even then often a difficult one. It is best approached by means of the 
reciprocal lattice, a concept which is not used in this chapter at all. A further section to the 
chapter on fibre photographs is certainly required. 
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The numerous diagrams are almost all excellent. Two or three could have minor 
improvements to make them more readily understandable. 

The book is very well produced and very few misprints have crept into the text. 

The index, while not very detailed, is sufficient, if one is prepared to look for a rather 
general heading and search a few pages of text to find the detailed reference required. This. 
is comparatively simple since the text is divided into small sections which have clear headings. 

The points of criticism raised above are comparatively minor ones which are almost 
inevitable in the first edition of a pioneering book. Future generations of students and 
teachers of x-ray crystallography will have reason to be grateful to the authors for the care 
and thoroughness with which they have made available their unrivalled experience in the 
wide field covered by this book. J. W. JEFFERY. 


Recommendations for Waste Disposals of Phosphorus-32 and Iodine-131 for 
Medical Users, by U.S. Department of Commerce, National Bureau of 
Standards. Pp.iv+11. (Washington, D.C.: Government Printing 
Office.) 
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ABSTRACTS FOR SECTION B 


The Texture of Diamonds used for Counting «, 8 or y Particles, as found from 
Divergent-Beam X-Ray Photographs, by H. J. GRENVILLE-WELLS. 


ABSTRACT. Since the radiation-counting efficiency of diamonds is a structure-sensitive 
property, X-ray experiments have been undertaken to see whether it is related to crystal 
texture. A technique for obtaining divergent beam photographs superimposed on Laue 
photographs has been developed, and results are given for 28 out of 38 colourless octahedral 
diamonds (edge 1-3 mm) whose ultra-violet transparency, $-particle counting efficiency and 
fluorescence are known. The relative strength of the anomalous extra x-ray reflections 
generally found for comparatively perfect diamonds has also been determined for these 38 
specimens. Counting efficiency does not seem to be closely related to texture and, in 
particular, mosaicity is not a necessary condition. 


Non-Linear Amplification in E.M.I, Photomultipliers, by J. F. RarF.e and E. J. 
ROBBINS. 


ABSTRACT. The output pulse heights obtained from E.M.I. photomultipliers are found 
not always to be proportional to the quantity of light incident on the photo-surface of the 
photomultiplier. ‘The photomultiplier ceases to be a linear amplifier of the incident light 
pulses when the charge density in the multiplier becomes too great. 

The relation between pulse height and energy loss of alpha-particles in various 
fluorescent materials has been determined at low values of charge density. 


On the Measurement of Current Density at the Cathode of a Glow Discharge 


and Aston’s Law at Low Pressures and High Discharge Voltages, by 
R. M. Cuaupuri and M. A, Baaqut. 


ABSTRACT. Amethod has been described for determining the current density on a plane 
cathode of a glow discharge over a wide range of pressure and voltage. It is shown that 
Aston’s relation, V=E+-F4/i/P, holds good down to as low a pressure as 2 x 10-? mm Hg 
and up to as high a voltage as 3000 volts. The work of the earlier investigators extends up 
to a voltage of about 1000 volts and a pressure not much below 1074 mm Hg. It is also 
seen from these experiments that the distribution of current over a plane cathode is a 
function of the discharge pressure and is in general independent of voltage and current. 


The Theory of Gaseous Arcs: I—The Fundamental Relations for the Positive 
Columns, by K. S. W. CHAMPION. 


ABSTRACT. It is postulated that the five fundamental relations, stated by Tonks and 
Langmuir to be sufficient to describe fully a low pressure positive column, also apply to high 
pressure positive columns. Previously, only two of these relations had been applied to high 
pressure columns. A ‘natural arc radius’ is postulated for low pressure arcs when the 
boundaries are removed an indefinite distance. Previous theories of low pressure arcs 
always assumed that the arc completely filled the cross section of its containing vessel and 
thus did not apply to low pressure arcs in large chambers or arcs at intermediate pressures. 

The theory of low pressure arcs with high electron densities in which ionization by stages 
is important is developed. It is then shown that by taking recombination into account this 
theory may be applied to intermediate and high pressure arcs with a high degree of ionization. 
A first approximation is also obtained for the theory of high pressure arcs in which thete is a 
low degree of ionization. All the theory presented includes the effect of a longitudinal 
magnetic field, and the no-field conditions are included as a special case. It is found that the 
effects of a magnetic field predicted by the theory are in agreement with experimental 
observations. 
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The Theory of Gaseous Arcs: II]—The Energy Balance Equation for the Positive 
Columns, by K. S. W. CuHaMPION. 


ABSTRACT. The theory of arcs is developed further by investigating the energy balance 
equation for high pressure positive columns. The energy balance equation put forward 
by Suits and based on the Nusselt relation is considered and some refinements in his analysis 
result in better agreement between theory and experiment. This theory only takes into 
account energy losses due to convection and conduction, but the neglect of radiation losses is 
justified for a wide range of conditions by reference to experimental values of radiation 
efficiencies. 

The effect of a longitudinal magnetic field on high pressure arcs is then investigated. 
From the point of view of the energy balance equation the action of a magnetic field is to 
reduce the contributions of the charged particles to the conduction and convection losses 
from the column. Horizontal and vertical arcs are considered separately, as the nature of 
the convection, and hence the effect of a magnetic field, is not the same in the two cases. 

The theory predicts that the action of a magnetic field will tend to a limit with sufficiently 
strong fields, and certain simple relations are deduced which should apply in this case. 
According to the theory, when the radius of a high pressure arc is reduced by a magnetic 
field the arc temperature must rise (incidentally, only by a relatively small amount) and the 
positive column gradient may increase or decrease, depending on certain relations which are 
specified. There are virtually no experimental data available for these conditions, but the 
scant data available are in general agreement with the theory. 


Observations on Radio-Frequency Oscillations in Low-Pressure Electrical Discharges, 
by N. R. Lasrum and E. K. Bice. 


ABSTRACT. The generation of high-frequency energy by electrical discharges in air at 
low pressures has been investigated by observations with several types of discharge tubes, 
at frequencies near 200 Mc/s and below 2 Mc/s. 

Random noise of high intensity was observed in both these frequency ranges, and was 
found to be associated with a field reversal in the discharge. 

Coherent oscillations were also observed at the lower frequencies, and there was evidence 
that these were the cause of the striations of the positive column. 


An Improved Frequency Equation for Contour Modes of Square Plates of Anisotropic 
Material, by R. BECHMANN. 


ABSTRACT. An improved approximation, based on previous work, has been derived 
for the frequency equation covering the four contour modes of vibration of square plates of 
anisotropic material. ‘The modes consist of three longitudinal and one shear mode. An 
X-cut quartz, a Y-cut EDT and Y%-cut lithium sulphate plates have been used as examples. 
The equation has been used to determine the elastic coefficients s,, and s,; of EDT and 
lithium sulphate monohydrate. 


Elastic and Piezoelectric Coefficients of Lithium Sulphate Monohydrate, by 


R. BECHMANN. 


ABS TRAC die A new determination has been made of the elastic and piezoelectric 
coefficients of lithium sulphate monohydrate and of their temperature coefficients. 


The Absorption Spectra of Single Crystals of Lead Sulphide, Selenide and Telluride, 
by A. F. Gipson. 


ABS TRAC T. The absorption spectra of single crystals of PbS, PbSe and PbTe have been 
examined in the temperature range from 20°K to 600°K. The spectra are characterized 
by a sharp absorption edge which coincides with the long wavelength limit of photocon- 
ductivity in each material. At longer wavelengths the absorption coefficient increases 
slowly with increasing wavelength, as in germanium and silicon. 


